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From the Editors 


A great deal of attention in the Soviet Union is being paid 
to elaborating on the scientific legacy of the founders of rocket 
and space science and technology. Since 1966 annual scientific 
lectures have been held in Kaluga, dedicated to elaborating the 
creative legacy and developing the ideas of K. E. Tsiolkovskiy ; 
since 1970 the F. A. Tsander memorial scientific lectures have 
been held in Moscow, Leningrad, Riga, and Dnepropetrovsk; since 
1971 various scientific and social organizations, the USSR 
DOSAAF [Voluntary Society for Collaboration with the Army, 

Air Force, and Navy] , and institutions of higher learning have 
been conducting lectures dedicated to some pioneers in the 
conquest of space. 

In 1977 the Presidium of the USSR Academy of Sciences 
adopted a resolution concerning the holding of a scientific 
session dedicated to the 70th anniversary of the birthday of 
S. P. Korolev. Participating in the work of this session of 
lectures were prominent Soviet scientists who worked with 
S. P. Korolev, researchers into the scientific creativity of this 
remarkable scientist, veterans of Soviet rocket technology, 
and the scientific societies of Moscow and other cities in 
the nation. 

The second set of S . P . Korolev memorial lectures took 
place between 29 March and 4 April 1978. About 50 reports 
were given at the plenary session and at the sessions of four 
sections. In addition, the symposium "The Role of the Researcher's 
Personality in the Genesis of Cosmonautics Theory and Liquid- 
Fuel Rocket Technology" was held within the framework of the 
symposium. 

Later, in connection with the creation of the USSR Academy 
of Sciences Commission on Elaboration of the Legacy of Pioneers 
in the Conquest of Outer Space (Corresponding Member of the 
USSR Academy of Sciences V. V. Raushenbakh, Commission chairman) , 
entrusted with study of the scientific creativity of Academician 
S. P. Korolev, Academician M. K. Yangel', Corresponding Member 
of the USSR Academy of Sciences G. N. Babakin, and other founders 
of space rocket science and technology, a resolution was adopted 
to conduct joint scientific lectures dedicated to the memory of 
the outstanding Soviet scientists and pioneers in the conquest 
of outer space. 

The 1979 scientific lectures were held in Moscow from 26 
February to 2 March. Over 1000 people took part in them, including 
prominent Soviet scientists and experts in space rocket technology, 
aviation and cosmonautics historians, veterans of Soviet rocket 
building, pilot-cosmonauts of the USSR, and talented junior 
scientists . 
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About 80 reports were given at the plenary and section ses- 
sions. They dealt with development of the scientific legacy of 
the founders of space rocket technology, problems of applied 
celestial mechanics and guidance, and topics in the theory and 
construction of aircraft engines and the history of aviation and 
cosmonautics. Three thematic sessions dedicated to the scien- 
tific creativity of A. M. Isayev, B. S. Stechkin, and S. A. 
Chaplygin were also held. 

The collection before the reader includes reports read at 
the sessions of the "Theory and Construction of Aircraft 
Engines" section and three reports on this subject given at the 
1978 lectures. The basic content of the collection consists of 
results of research on various aspects of the theory of opera- 
tion of jet engines at high supersonic speeds: high-temperature 
gas-turbine jet engines (GTD) , hypersonic ramjet engines (GPVRD) , 
and liquid-fuel rocket engines. 

Discussions involve the problems of testing modern complex 
high- temperature GTD and trends in the development of testing 
methods, and requirements for the principles of testing future 
aircraft engines are formulated. Methods of selecting the 
basic parameters of cooled high-temperature turbines of aircraft 
GTD are examined. 

An approximate model of the operation of high-speed PVRD 
[ramjet engines] (SPVRD [supersonic ramjets] and GPVRD) is 
examined, the results from analysis of errors in experimental 
determination of these engines' parameters and results from 
computational study on selecting the optimal geometry of the 
GPVRD air inlet are cited, the effect of irregularity of a 
supersonic flow with shock waves on friction and heat exchange 
in the GPVRD duct is investigated, and the results of numerical 
calculation of plane supersonic flows in ducts are cited. 

An analysis is made of means of hydrodynamic modeling of 
operation in ZhRD [liquid-fuel rocket engines] in the launch 
mode; the hierarchy of mathematical models of aircraft control 
ZhRD is examined. 

The collection is aimed at the use of specialists working 
in the field of reaction engines for aviation and space systems. 
It may also prove useful to scientists, teachers, and students 
at aviation institutes . 
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THEORY AND DESIGN OF FLIGHT-VEHICLE ENGINES 
(Proceedings of the Joint Scientific Lectures 
Dedicated to the Memory of Soviet Scientist-Pioneers in the 
Conquest of Outer Space: Moscow, 26 February-2 March 1979) 

V. T. Zdanov and R. I. Kurziner, Eds. /5* 

CURRENT PROBLEMS IN TESTING AVIATION ENGINES 1 

V. M. Akimov 


The new and greater capabilities opening up before designers 
and researchers because of the introduction of computers and 
assorted powerful numerical methods are in many cases creating 
the illusion that the role of experiment is diminishing, especially 
in the design of aviation and space technology. For aircraft 
engines the problem of testing is quite an acute one and requires 
answer© to many new and complicated questions. We shall examine 
this subject using the example of aircraft jet engines, bearing 
in mind that the problem is characteristic of rocket engines 
as well — with the difference, of course, produced by the 
specific features of both engine types. 

Gas-turbine jet engines began being produced when the theory 
and calculation methods of engine and element operation had 
already been basically developed (in our country this was accom- 
plished primarily by Academician B. S. Stechkin and such 
other scientists as Prof. B. V. Uvarov, Prof. K. V. Kholshevnikov, 
and their scientists) . Fairly refined methods of calculating 
gas-turbine engine strength also developed rapidly. Despite 
this, experimental work on operational development and testing /6 
of engines under construction did and continues to constitute 
the basic content of gas-turbine engine construction with 
regard to time and cost. One of the major reasons for this 
situation is the inadmissibility of any markedly excess weight 
and size for engines and aircraft components in general. This 
necessitates experimental testing of the adequacy of selected 
performance margins for engine components and systems, which 
in most cases are notable for high load levels. 


Numbers in the margin indicate pagination in the foreign text. 

1 In preparing this report the author employed articles from various 
authors published in Ispytaniye aviatsionnykh dvigateley [Aviation 
Engine Testing], UAI Press, Ufa: No. 1, 1972; No. 2, 1974, No. 3, 
1975; No. 4, 1976; No. 5, 1977; No. 6, 1978. 
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Data characterizing the process of experimental operational 
development of aviation GTD [gas-turbine engines] are given in 
Figs. 1-3. The many years' experience of engine-building firms 
indicates that the period of operational development tests 
averages 50-60 months (Fig. 1) . Note that 30 to 50 months pass 
from the first test-bed tests to the first test flight. Several 
dozen experimental engines are used in the operational develop- 
ment process. A look at Fig. 2 reveals that there is a tendency 
towards reduction of their number , but the total accrued opera- 
ting time is not decreased (Fig. 3) because the average accrued 
operating time for each example is increased. Out of 10^ hours 
of total accrued operating time, (1 + 2) *10^ hours are spent on 
engine operation under simulated field conditions in altitude 
chambers and wind tunnels. The total cost of operational devel- 
opment for an experimental engine is 130-150 times greater than 
the cost of full-scale production of one engine model in 
sustained production. This ratio is 30/60 for the airframe. 
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Fig. 1. Experimental operational development periods 
for various types of GTD in months (before beginning 
operational use) from the year operational development 
begins . 


— trend line. 


Development of numerical methods based on use of powerful 
modern computers led to an abrupt increase in assorted para- 
meters, defined in computed design, characterizing operation of 
all GTD components. The total number of computation operations 
in the design of an engine has increased from 10^ to 10-*-^ in 
30 years. Universal methods allowing high-precision computation 
of the stressed state of extremely complicated parts, methods 
of optimizing an engine and its components at all levels of 
design (engine — in the aircraft, assembly -- in the engine, 
part — in the assembly), simulation models, such as operations 
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Fig. 2. Number of experi- 
mental GTD models relative 
to year operational develop^ 
ment begins. 


and necessitated involvement 
for its solution. 


models, allowing, in particular, 
tracing the possible consequences 
of various failures, and other 
such possibilities presumably 
should have simplified the tasks 
of tests and reduced their required 
number. This did not occur in /8 
practice, however. The complica- 
tion of engines and their require- 
ments, occurring together with 
refinement of computation and 
design methods, was so significant 
that experimental operational 
testing tasks also became constantly 
more complex. In the last few 
years a number of sometimes con- 
flicting factors have made the prob- 
lem of engind testing rather acute 
of new approaches and new methods 
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We shall examine the 
factors determining the cur- 
rent requirements for avia- 
tion engine tests; we shall 
deal first with the factors 
requiring increased volume 
and complexity of data 
obtained in testing. 
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First among such fac- 
tors is the constantly 
increasing design and para- 
metric complication of GTD. 

We shall use simple examples 
•♦STv.years to show the effect of com- 
plication of engines on 
the complication of their 
tests. For instance, in 
modern GTD the temperature 
level of the gas ahead of 
the turbine is approximately 
1600K (it will reach 1800- 

1900K in the next generation of engines) , so only cooled turbine 
blades, notable for construction complexity, are employed. A 
cooled-blade thickness variation within 0.1 mm leads to a 
wall-temperature difference of v50K, impermissible with relatively 
small strength margins from the standpoint of required blade 
reliability and lifetime. Since such technological variability 
cannot be precluded in experimental models, the wall temperature 
of each operating blade has to be measured during operational 


Fig. 3. Total running time of GTD 
during operational development 
relative to year operational devel- 
opment begins. 
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testing. On the other hand, all modern GTD are two^ and three- 
shaft. Obviously, it is practically impossible to tap sensor 
signals from the internal rotor using traditional current- 
collectors. Consequently, fundamentally new means of gathering 
data are needed. Furthermore, the specific unit weight of modern 
engines has reached a level on the order of y sw = 0.125, corres- 
ponding to a thrust/weight ratio of 8. This requires designating 
minimum performance margins, and thus maximum experimental testing 
of the construction. In sum, we have to obtain a constantly 
increasing volume of information, which at the same time becomes 
more and more difficult to obtain. 

Second among such factors is the complication of use condi- 
tions, including engine operation for most of a flight with inlet 
air parameters above atmospheric on the test-bed (an SPS [SST] 
engine) and integration of the powerplant components and airframe/9 
(for example, the powerplant of V/STOL [vertical/short take-off 
and landing] aircraft with a surface-blown wing) . Both these 
factors require testing to provide complex simulation of flight 
conditions on special, expensive and power-consuming stands. 

Third among the factors complicating GTD testing would be 
the extremely high reliability of modern engines. For passenger- 
air craTF-engines^ firmly established in operation and production, 
this is characterized by mean times-between-f ailures-in-f light 
of (3 + 5) *10^ hours, and sometimes even more, and for some 
specific failures, by values on the order of 0.5*10^ hours. Con- 
firmation of this reliability at the testing stage using tradi- 
tional statistical methods with any acceptable trustworthiness 
would require a multiple increase in engine running time during 
operational testing, and a very large increase, as seen in Fig. 3. 
The practical impossibility of such a formulation of the question 
led to the need to develop and introduce special indirect methods 
of evaluating engine-component and system reliability, before 
launching full-scale production, in a relatively short time of 
testing . 

Fourthly, difficulties arise in life testing because of the 
very large lifetime values of many modern engines. For passenger- 
aircraft engines the specified life (added up from several over- 
haul periods) is 3-l()4 hours, and this is now common. If the 
engine performance for the entire specified life were checked 
on a full time scale, operating ^10 hours a work day, such life 
tests would take about 10 years. The senselessness of this 
situation made it necessary as engine lifetimes increased to 
develop and employ new, accelerated methods of testing engines, 
so-called equivalent tests. 

Fifthly, new requirements for engines, associated with 
environmental protection, led to the need for GTD tests new to 
the engine-builder — tests to determine the content of harmful 
materials in the combustion products exiting the engine and to 
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evaluate engine noise. New testing methods and new equipment 
were needed. 

Sixthly, the implementation of "Standards for Airworthiness 
of USSR Civil Aviation Aircraft" (NLG SSSR [USSR Airworthiness 
Standards] ) , which met the requirements of the International 
Civil Aviation Organization (ICAO) , necessitated some review of 
testing procedures and methods testing the characteristics of /10 
new engines. Some new types of tests (tests on autorotation 
without lubrication, for instance) were introduced as part of 
these tests, the results of which were used to assess the suita- 
bility of an engine for flight operations, and the programs and 
procedures of certain tests had to be reviewed (e.g., the two- 
hundred-hour tests, traditional for our country, were replaced by 
one-hundred-fifty-hour tests, which better meet the generally 
accepted international requirements) . 

Along with the above factors necessitating expansion of 
experimental data, tendencies are found to hinder meeting of these 
requirements. First of these is the constant increase in the 
cost of experimental developments, associated with the compli- 
cations of engines. This cost has increased ^20 times since the 
first gas-turbine engines were built. Obviously, this circum- 
stance forces us to try to reduce the number of experimental 
models used up in the experimental operational development of 
an engine. 

Secondly, requirements for fuel and energy-resource economy 
face us with the task of constantly finding possibilities of 
decreasing the extent of testing (including in special high- 
altitude installations) . 

Thus, two conflicting groups of factors are making it neces- 
sary to obtain increasingly extensive data on a decreasing number 
of experimental models and in a shorter time of testing. The 
practice of engine-building has now developed certain approaches 
to solution of these problems, considering the above conflicts. 

A number of trends have been mapped out in order to provide 
effective solution of the engine-testing problem. Here are 
some of them: 

— introduction of non-contact measurements and combining 
of tests; 

— planning of the experiment; 

— extensive automation of the experiment and use of mathe- 
matical models of engines in tests; 

— production of special test-beds simulating an increasingly 
wide range of flight conditions; 
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— advancement of accelerated equivalent and cyclic tests 
for testing lifetime? 

— rejection of certain types of testing; 

— development and use of special certification tests for 
engine operational reliability; 

— development and introduction of new types of tests (for /II 

noise and emissions) . ’ 

We shall deal briefly with the above topics, and we shall 
go into somewhat greater detail on some which, in our view, may 
be of some interest also beyond the field of aviation gas-turbine 
construction . 


Non-contact measurements. Various methods of non-contact 
measurement of GTD parameters have become widespread in recent 
years. We shall examine several such methods. 


One interesting instrument for non-contact measurements is 
that which has been termed ELURA (cathode-ray device for recording 
blade- vibration amplitudes). The instrument's operation is based 
on the discrete-phase method. A fixed pulse-transducer (2 in 
Fig. 4) generates a pulse at the time some stimulus with an 
angle pitch equal to the blade's passes by it. Such stimuli may 
be holes, teeth, etc. Another pulse-transducer (1 in Fig. 4) 
generates a pulse when the tip of a blade whose vibrations are 
being measured passes by it. The time interval between pulses 
of the actuating and synchronizing transducers Ti m p is proportional 
to the path transition time hi. (see Fig. 4) 


imp 


= /u., 


(D 


where u , is the circumferential velocity at the rotor periphery. 
In turn, linear marks made by an electron beam on the screen 


of the recording device hi '. are proportional to x 


hi ! v hi. 
x % 


( 2 ) . 


imp' 


so 


If there are no blade vibrations, the hli degenerate into 
points. Several new types of instruments have been created with the 
use and in the development of this principle. Both transducers 
can be placed in the housing above the rotor. By placing the 
transducers over different points on the blade face, measurements 
can be made of different forms of vibrations, including some 
torsional ones. 


Use of ELURA-type instruments may not always be effective, 
however, e.g., when studying high-frequency blade vibrations by 
certain complex shapes. In such cases it remains necessary to 
read signals from transducers mounted on the internal rotors of 
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Fig. 4. Schematic diagram of blade-vibration 
amplitude measurement using the ELURA instrument. 

1 — actuating pulse transducer, 2 — reference 
synchronizing pulse transducer, 3 — cathode- 
ray recording device. 


multishaft GTD. Then we have to resort to microtelemetry , i.e., 
use of miniature transducers that rotate together with the rotor. /13 
Requirements for them are very strict: for example, one adopted 
model provides for operation at a centrifugal acceleration of 
3*1()4 g and allows a transducer temperature of 125 °C. Use of 
microtelemetry is one of the developing trends in non-contact 
retrieval of data during engine tests. 

Optical pyrometries (Fig. 5) have recently become widespread 
for measurement of operating blade temperatures. The operating 
principle of the optical pyrometer is quite familiar. Difficulties 
encountered when using this instrument for aviation GTD measure- 
ments are also fairly characteristic. In particular, they are 
associated with the presence of soot in the flame. In most cases 
the error (when measuring blade-wall temperatures up to ^1000°C) 
can be considered not to exceed approximately 10°C. Use of a 
scanner allows measuring the temperature of each operating blade, 
something practically impossible when measuring with thermocouples. 
The lower right corner of Fig. 5 shows an example of the type 
of recorded data: the temperature of each of z blades is measured 
separately. 

Note that research is under way on use of optical methods 
such as Raman spectrography , which should allow, in particular, 
measurement of the temperature of any level at a point in a 
stream of high-speed gas. 

We shall mention two other promising methods of non-contact 
measurement. First of these is vibroacoustic analysis, including 
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Fig. 5. Non-contact measurement of operating 
turbine-blade temperatures using an optical 
pyrometer . 


1 -- blade, 2 — turbine housing, 3 — objective 
lens, 4 — optical fiber, 5 -- photodetector, 

6 -- amplifier, 7 — recorder. 


analysis of noise emitted by the engine. In this case, micro- 
phones mounted outside the engine are used to draw an acoustic 
portrait in different operating modes. Noise analysis can, in 
principle, be used to evaluate, for instance, closeness to surge, 
bearing defects, breakage of a small piece of blade, and many 
other data. One difficulty is separating the very weak signal 
when there is a very high general noise level. Analysis of the 
noise spectrum requires complicated and sometimes narrowly 
specialized equipment. 

Also being developed are methods of laser Doppler velocity 
measurement (LDIS) , which allow detailed study of flow structures. 
Some successes have already been achieved in using LDIS to study 
flows, even those of complex structure, such as the flow in the 
channels of an axial-flow compressor. 

Non-contact measurements not only provide information diffi- 
cult to obtain with other methods, but also permit combining /15 
GTD tests with different objectives. For example, tensometric 
measurement of internal-stage blades, transmitting sensor signals 
using microtelemetry, study of low-pressure stage blade vibra- 
tions using the ELURA instrument, pyrometer measurement of 
turbine blade temperatures, reading of the engine's acoustic 
portrait, etc., can all be accomplished simultaneously. Conse- 
quently, non-contact measurements are a powerful method for 
meeting the conflicting requirements placed on testing of modern 
avaiation gas-turbine engines . 
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Experiment planning. Experiment planning methods recently 
finding use in aviation engine-building (previously employed in 
other branches of technology) have provided an appreciable effect 
in solution of problems arising in testing modern GTD. The 
theory of experiment planning is, of course, based on the concept 
of a multifactor experiment, while the traditional experiment 
in tests of engines and their components was always a single- 
factor one. Combination of factors (a?i, X2, ...» in measure- 

ments is chosen in accordance with certain rules, and processing 
of the experimental data produces regression equations of the 
y = f (x i, X2 , . .., a?k) type, which associate the characteristic 

being sought with the factors determining it. 

Experiment planning allows, first of all, reducing the 
required volume of experimental data with a given accuracy of 
results and, secondly, it opens up great possibilities for 
optimizing the system being tested, e.g., for finding the most 
suitable adjustment of engine components. Figure 6 illustrates 
the former using as an example the results of expensive tests 
to read the high-altitude characteristics of an engine in an 
altitude chamber. The graphs in this figure show that ten points 
were needed to draw the throttle characteristics at five alti- 
tudes (instead of the minimum thirty with a single-factor 
experiment) . The agreement of characteristics obtained using 
regression equation R = /(H, n) , where H is the flight altitude 
and n is the r.p.m. , with the results of control single-factor 
experiments for two values of H is good. 

The idea of one method of optimization (simplex) is shown 
for a two-factor (a?i and a^) experiment in Fig. 7. First, 
experiments corresponding to combinations of factors at points /17 
1, 2, and 3 of the factor plane are set up. The vertex of the 
simplex (point 2) where the worst value of parameter y was 
obtained, is mirror-reflected and coordinates and a?2 of new 
vertex 4, for which the next experiment is performed, are found. 

The points of the new simplex 1-3-4 are used to select the 
worst one according to the results, and it is again mirror- 
reflected (point 5), etc. In most cases the simplex method 
allows entering the optimum region in a few steps. Most effec- 
tive are such optimization methods in experiments whose results 
had previously been sometimes difficult to assess even roughly, 
in the operational development of combustion chambers, for 
example . 

Experiment automation and use of mathematical models of 
engines . A great deal has been written about experiment auto- 
mation, in connection with which we shall only briefly mention 
one aspect of this question. We note that when discussing auto- 
mation we most often mean automation of the collection and 
processing of measurement results, which earlier had been per- 
formed manually. Labor productivity thus increases and the 
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Fig. 6. Throttle characteristics 
of an engine at M = const, and 
H = var. , obtained using regres- 
sion equation ( ) . 

o — combinations of factor values 
used to conduct experiment. 
o> — experimental points obtained 
in control single-value experiments. 



Fig. 7. Diagram of 
optimization with 
simplex method. 


efficiency of equipment use 
is enhanced; the time con- 
sumed by the experimenter 
in making decisions is 
reduced if processing 
takes place at the rate 
of the experiment. There is, 
however, no sharp reduction 
in the volume of testing 
if the approach to organi- 
zation of the experiment is 
not altered. A decisive 
effect can be obtained 
with a fundamental change 
in methods of obtaining data 
about the item being tested 
and in methods of evaluating 
: the test results. For ex- 

racter iqtics ample, if in a multi-level 

. , automation system powerful, 

.j nrr r(lftrM _ high-level computers compare 

g g the experimental data with 

various mathematical models 
factor values (thermal gas-dynamic, dy- 

iment. namic, etc.) of the test 

nts obtained engine, fairly complete 

ue experiments. identification of the 

engine with its models can 
be accomplished using the 
results of a limited number of measurements. 
Use of such methods of evaluating test 
results according to limited information, 
when models of an engine or one of its 
assemblies are matched with experimental 
material, allows reducing the required vol- 
ume of experimental data. Prospectively, 
with use of non-contact measurements, com- 
bination of goals of different tests, and 
use of a set of mathematical models in an 
automated system we can expect a marked 
reduction in the volume of tests done to 
assess the achieved parameters of the 
assemblies, and of the entire engine, and /18 
a choice of directions to take in further 
operational development. 


Simulation of flight conditions. The expansion of engine- 
use conditions and the increase in air-mass flows demand con- 
stant complication and increased cost of special installations 
for altitude-velocity tests. Thus, one of the main directions 
to take in bringing these requirements into line with economic 



ones is employment, in all possible cases, of simplified experi- 
mental installations for partial simulation of flight conditions. 
We have in mind, for instance, simulation of velocity at slowed 
air parameters, heating of incoming air by mixing it with the 
exhaust gases of another engine operating on an adjacent test- 
bed, use of pulsed tunnels, etc., etc. This is a huge area, 
fairly specific for aviation engines, and there is no reason 
to dwell upon it in more detail within the framework of this 
article . 

Equivalent and cyclic tests. Here the main question is 
that of equivalence (with respect to failure rate) of equivalent 
test conditions to operational conditions. We shall examine 
some possible approaches to assessing equivalence of GTD 
conditions . 

1) Often it is expedient to employ the measure of lifetime, 
associated with strength margins. Then the principle of equi- 
valence is reduced to equality of safety margins in accelerated 
and operational tests . Since we can consider that the lifetime 
of a GTD is determined primarily by the turbine blades and disks, 
the equality of statistical safety margins at increased tempera- 
ture is adopted as the measure of the equivalence of conditions. 
At each temperature T for a given material we can write equation 

ff t *• C , (3) 

where t* is the prefailure time with exposure to static stress d . 
The equivalence principle, based on equality of safety margins, 
is illustrated in Fig. 8. If Ti, t i, and 6i are the temperature, 
operating time, and effective stress under the conditions simula- 
ted by the equivalent tests, and the long-term strength 

limit at these T. and t., then the safety margin in this condi- 
tion 

k; /*- (4) 

If in another (t - m) condition tj is the operating time 
with which the safety margin 



equals that in the ith condition, then this (1st) condition is /19 
equivalent to the ith. Using the statistical failure model (3) 
and a linear model of damage summation we can obtain the equi- 
valence condition 

^ea’^ * « (€) 

where t^ e q is the equivalent time in the 1st (selected as equi- 
valent) condition and n represents all conditions simulating 
the equivalent. Equation (6) provides a decisive curtailment 
of testing time. For example, if we double the operating time 
of an aviation engine in the maximum condition, in many cases we 
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Fig. 8. For evaluating the 
equivalence of long-term 
strength conditions. 
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can consider that, by vulnerability 
to damage of the hot portion, this 
is practically equivalent to a 
lifetime of engine operation. 

2) Widely prevalent are cyclic 
tests in which equality of the number 
of cycles of low-endurance fatigue 
in a lifetime and in cyclic tests 
is taken as a measure of equivalence. 
The effect of other conditions on 
failure rate of the hot area and 
rotor of a GTD is considered less 
substantial. Loading cycles of the 
cycle type shown in Fig. 9 for a 
passenger aircraft TRDD [fanjet] 
engine are reproduced in the tests. 
The short duration of the cycle 
allows reproducing many thousands 
>f them in a short time. Normally 
in cyclic tests the statisti- 
cal scatter of the results 
of low-cycle tests is taken 
into account. Based on the 
possible difference in the 
properties of materials, 
tests are done with a margin 
against the number of start 
and shutdown cycles per 
lifetime : 


N 


= cn /n 
res 


(7) 


Fig. 9. 
cycles . 


& io iz m t.min 


Examples of TRDD test 


1 -- idle, 2 — take-off mode, 

3 — cruise mode, 4 — descent 
mode, 5 — idle in flight, 

6 — thrust reversal, 7 — inter- 
mediate idle, 8 — shutdown. 


test 

where N reg is the number of 
cycles per lifetime, Nt es t/ 
per test; C is the margin 
by number of cycles consider- 
ing scatter, and n is the 
correspondence factor of the 
operating and test cycles. 


The value of C is deter- 
mined using the scatter field 
of the results of tests on 

prototypes (see Fig. 10). If a crack appears in a disk at 
N t est number of cycles, it is assumed that this is the best disk. 

The worst disk at the same load fails at N' cycles (see Fig. 10) . 
Hence, C = N-^est/N' . The value of coefficient C can in practice 
approach C = 2-3. Also possible is an overstress version of 
testing. If a bad disk fails at N' cycles and cl^ stresses 
(Fig. 10), the best runs N' at stress a"l- Then with N' operating /21 
cycles we can avoid increasing the number of test cycles (N. .j. =N'), 
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scatter 

field 


log hi' log A/ log V . 

Fig. 10. For accounting for 
scatter of results of low-cycle 
tests in GTD cyclic tests. 


but increase it C = o'\/o\ times 
in stress tests. 

3) The statistical approach 
to assessing equivalence of 
conditions reduces essentially 
to equality of probabilities 
of failure-free operation in 
equivalent tests and in a 
lifetime: 

(g^ 

where P(ti, zg) is the proba- 
bility of failure-free operation 
in time tj_, in mode Zj_, and 
P(t 2 , Z 2 ) is the same in time 
t 2 in mode z . 



indices, we write 
t 

where X C p = / X(t)dt 
where T is the mean 
instead of (9) 


Using the mean statistical 

A 9 (x,)t 4 , (9)! 

is the mean failure rate. Since X CD = I/T, 
running time per failure, we can write 

t,/T(x„) s t t /t (i t b (10) 


Studies by the author and A. A. Morozov found statistical 
relations linking engine running time per failure to its opera- 
ting conditions, especially flight durations ^tf^) and propor- 
tions of heavy duty conditions per lifetime (z) . Specifically, 
for a large number of engines it was shown that 

a 

t - t fl a . (II) 

Substituting (10) in (9), we find that if subscript 1 in (9) 
means equivalent tests and subscript 2 means flight: 


^ec / V s ( ^test/^fl ^ (*test/'^ (j2) 

where x is the engine lifetime, is the duration of the 

test cycle, and ztest is the relative proportion of heavy duty 
conditions in the tests. From (12) it follows that the greater 
the proportion of heavy duty conditions and the shorter the test 
cycle, the shorter the time of equivalent tests. 


Note that, ignoring the second relation in the right part 
of (12) , bearing in mind the small value of the exponent in it 
and considering the closeness to unity of the exponent in 
ratio t t /t fi , we can write as an approximation 
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^ eq Atest? 1 ^ /tf 1 . • 


(13) 


But t e „/t t ^ st - N test , and i/tf^ = Nf^. Then from (13) we obtain /22 


the condition for equivalence of conditions 

N, . = N-, 
test fl 


(14) 


i.e., we come to the idea of cyclic testing. 


Practical use of the statistical approach to equivalent tests 
is difficult because of the absence of models of the (11) type 
for new engines. Compromise approaches, combining equivalent 
(by long-term strength of combustion area) and cyclic tests, are 
usually used in practice. 


Rejection of some tests. We always have to seek out possi- 
bilities for reducing the volume of tests and to review the 
range of GTD tests. For example, in many cases it has proven 
possible to eliminate some tests in mass-production plants, re- 
placing them with some kind of control. Commission tests, as a 
means of monitoring control over the stability of production of 
certain engines long in production, can be replaced by other 
methods of stability control, e.g., dynamic vibration tests of 
blades, etc. Two-step routine delivery tests using modern methods 
of diagnosing engine condition are replaced by single-step ones 
in many cases. Without dwelling further upon these questions, 
we note merely that, in addition to the necessary complication 
of the testing system, we must constantly work to make it smaller 
and simplify it in all cases where this is possible. 


Special tests to check the operational reliability of engines. 
The system of tests to check indirectly the reliability of an 
engine and all its systems has been developed, is in operation, 
and is constantly being improved. The requirements for these 
tests have been formulated in the airworthiness standards, and 
their positive results are a mandatory condition for engine cer- 
tification . 


The basic idea of tests is that they are a brief check of 
engine components and systems under assorted influences and 
operating conditions in order to determine the presence and 
adequacy of performance margins. 

Three questions always arise when developing and refining 
such a system of tests: what to test (test list), how to test 
(methods and programs) , and how to assess the results of the tests. 

The main, but not the only condition for formulating the /23 
list of tests is checking of performance margins as applied 
to typical failures (for example, fatigue failures of blades , 
surges, icing, foreign-object damage to the compressor, etc., etc.) 
or rare but more dangerous ones (disk failure, for example). 
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The system of testing now in effect contains several dozen types of 
tests. Therefore, despite the relatively short duration of most 
tests, a full set of them takes one or two years and requires use 
of several engines. Some tests, such as determination of blade 
vibration characteristics, checking of gas-dynamic stability 
margins, and examination of the lubrication system, blend into 
a whole set of operations including several engine overhauls, 
complex preparation, and many placements of the engine on test- 
beds, including in altitude chambers. 

The decisive question when developing test procedures is the 
choice of the level of loads (conditions and external influences) 
that the engine should withstand and at which its performance 
margins should be checked. For example, with strain measurement 
of turbine blades (operating) and compressors (operating and 
turning) we determine the effect on vibration stress level of 
such factors as greatest possible irregularities of total air 
pressures in operation, maximum possible working medium pressures 
and temperatures under high-altitude and -speed conditions, etc. 

The strain measurement should be performed throughout the entire 
range of r.p.m.'s, from idle to maximum, since both resonance and 
stall vibrations may arise at any r.p.m. values, as shown in 
Fig. 11, for example. 

Testing, for instance, 
of anti-icing system effective- 
ness should, in accordance with 
adopted standards, be performed 
under air- temperature and arti- 
ficial icing conditions cor- 
responding to the graph in 
Fig. 12. In nature it is very 
rare, but possible, to encoun- 
ter icing conditions more 
severe than those shown in 
Fig. 12. It is suggested that 
when such rare situations do 
develop in some areas, flights 
into them should be avoided 
at that time, since protection against icing is extremely diffi- 
cult under these conditions . 

We shall examine one other case: testing of engine performance 
after a bird strike. These tests check engine operation when sev- /25 
eral birds of average size (weighing 0.5-0. 7 kg) and small size 
(weighing up to 115 g) or one large bird (weighing at least 1.8 
kg) strike the engine. These conditions were selected based on 
analysis of worldwide experience in operating aircraft with GTD. 

The tests are conducted on an installation whose diagram is shown 
in Fig. 13. An air-driven gun has to be used to impart to the 
bird enough speed so that its velocity of impact with the engine 



Fig. 11. Sample variable stress 
curve . 


21 



equals the aircraft's cruising 
speed at the altitude where 
a bird strike is still possible 



Fig. 12. Recommended curves of 
water content for testing anti- 
icing systems. 

1 — prolonged action, 2 — 
action lasting up to ^1 min. 


Assessing test results 


requires a very well-founded 
regard for the consequences of 
some influence (if it takes 
place) and for the level of 
performance margins established 
by the tests. True, in some 
cases the assessment of conse- 
quences is unambiguous and 
causes no difficulties — and 
in these cases the assessment 
of test results proceeds accor- 
ding to the "yes or no" prin- 
ciple. Such is the situation 
in tests on housing impenetra- 
bility during blade breaks — 
if the broken blade does not 
pierce the housing and all 
secondary damage is localized 
within the engine, the test 
results are positive. 


3 c An opposite example 

can be served by assess- 

^ ment of damage from 

. - > a bird sucked into the 

engine. Some damage is 

unavoidable under these 

circumstances, but which 

to consider tolerable 

is not obvious without 

analyzing later flight 

, 0 ~ , i* . . . , situations. The follow- 

Fig. 13. Bird-launching installation. . . . , 

3 3 ing assessment criteria 

1 — GTD being tested, 2 — pneumatic have been adopted. A 

cannon, 3 — bird, 4 — compressed bird weighing 1.8 kg, 

air bottle. after striking the 

engine, should not inflict 

damage dangerous to a multi-engine aircraft. Spontaneous or 
forced shutdown of one engine is not viewed as dangerous to the 
aircraft. Impact of several medium-sized or small birds should 
not lead to engine shutdown or to inadmissible (for flight con- 
tinuation) deviations of parameters, since small birds may hit 
several engines. 




Tests measuring engine noise and emission of toxic materials. 
These types of tests, comparatively new for aircraft-engine build- 
ing, have given rise to a number of new questions. The major one 
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for noise measurement is that tests to determine GTD acoustic 
characteristics require creation of a free acoustic field, i.e., 
these tests are impossible on common test-beds. There should be 
no noise-reflecting structures by open test-beds for acoustic 
measurements, even 100 m from the engine. If service structures 
are present at distances up to 100 m from the line of the engine 
axis, their walls should be covered with sound-absorbing material./28 

In tests measuring toxic emissions the difficulties are 
associated with some equipment and procedural questions. First 
to be determined in these tests are engine exhaust concentrations 
of nitrogen oxides (NO x ) , carbon monoxide (CO) , unburned hydro- 
carbons (HC) , and smoke. Figure 14 shows examples of change in 
the concentrations of these compounds relative to engine settings. 
These concentration levels, in g per kg of fuel, indicate that 
quite infinitesimal contents of these compounds, relative to the 
total volume of engine exhaust gases, have to be measured. For 
example, the required range of NO x measurements is 0-0.1%, 
of HC, 0-0.2%. On the other hand, there are ICAO-recommended 
standards for total content of toxic substances emitted in an 
airport region over a rather prolonged standard cycle (Fig. 15) , 
including taxiing and waiting on a taxiway, take-off and climb 
to ^9 00 m, descent from this altitude, landing run, and taxi 
to park. Considering the changes in the concentration of hazardous 
materials with engine settings, this requires exceptionally high 
accuracies of instantaneous measurements of these concentrations. 
Modern methods have been developed to detect up to 10“ 6 to 10“7 
parts-by-volume of toxic substances in exhaust gases. There 
was a need, however, to employ different measurement principles 
for each compound: chemiluminescence for NO x , flame-ionization 
for HC, and a method based on nondispersed infrared absorption 
for CO. As a result, complex and expensive equipment has to be 
used in tests to measure emissions of hazardous substances. 

Engine complication has led to the fact that the total cycle 
for producing aviation GTD, including time studying separate 
components, developing the experimental prototype, and, finally, 
building and developing an experimental model put into large- 
scale production, has stretched to 15 years. And experimental 
work continues to play a dominant role at all stages of this cycle. 

* * * 

To sum up, we can formulate certain conclusions from the 
above material. 

Firstly, the complication of aviation engines and require- 
ments for them are outpacing progress in design and calculation /29 
methods, leading to complication and increase in the volume of 
tests. Thus, we need to have further intensive and purposeful 
improvemeht of, on one hand, design methods and, on the other, 
means and methods of experimental engine development. Otherwise 
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Fig. 14. Example of change in the concentration 
of toxic substances with engine setting. 


andin 


Fig. 15. Conventional standard "take-off - 
landing" cycle, for which total emission of 
toxic substances is determined. 

1 — taxiing, 2 -- take-off, 3 — climb, 

4 — descent. 







a further increase in conflicts between requirements and the 
capabilities of experimental operational development is un- 
avoidable. 

Secondly, in connection with the very large time period 
between the beginning and end of the engine-development cycle, 
there has to be very reliable prediction of a wide range of 
required engine properties, so as to reduce to a minimum the 
development of changes whose need develops and is discovered 
in the process of developing and testing an engine and its 
components . 

Thirdly, further, long-term development of engines merely 
through constant complication of them is impossible. An active 
search must be undertaken for fundamentally new solutions 
involving the work process, design-technology solutions, 
and materials. 
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ERRORS IN EXPERIMENTAL DETERMINATION OF 
HIGH-SPEED RAMJET PARAMETERS 


Yu. M. Shikhman, M. V. Strokin, V. A. Uvarov 


One of the pioneers in cosmonautics and rocket technology, /30 
F. A. Tsander, focused considerable attention on refining the 
thermodynamic cycles of jet engines and pointed out the promise 
for use of VRD [jet engines] in high-speed flight [1, 2] . Some 
confirmation of his prediction is found in ramjet engines designed 
for use as powerplants of aricraft flying at Mach numbers from 
3 to 6 and more [3-5, 7], 

Now, of course, serious attention is being paid to experi- 
mental research on the operation, parameters, and characteristics 
of GPVRD [hypersonic ramjet engines], and tests of GPVRD 
on flying laboratories are under discussion [5] . The complexity 
of conducting tests at high airflow M, both on ground-based 
test-beds and on flying laboratories, is what lends importance 
to the study of possible errors in experimental determination of 
high-speed ramjet parameters. The results of such research may 
be of help in choosing more accurate methods of processing experi- 
mental data and developing requirements for measurement accuracy. 

In this article we attempt to undertake such research. 

We use as an example a GPVRD operating at M 3 — a so-called 
dual-mode GPVRD (DGPVRD) [3] . One possible design for a DGPVRD 
flow duct is shown in Fig. 1. Its distinguishing feature is that /31 
the fuel combustion process can take place in a sub- or supersonic 
flow, depending on the fuel consumption or M. Accordingly, we 
shall call these the sub- and supersonic combustion modes. 

The basic parameters measured in tests of these engines are 
the static pressures Pf, measured on the walls of the flow duct, 
the temperature of these walls T w ^ , and data on fuel consumption 
G-p- Assume we also know the static pressure of the incident air 
flow Pjj, its stagnation temperature Tg, and the pressure behind 
a normal shock PA ' . In addition to these three parameters, we 
may also know others which also characterize the energy factors 
of the air flow, e.g., Pj$, Tg, and Pg' [sie] . 

We need the measured parameters to determine engine and 
incident-flow parameters: incident-flow Mach number Mh and angle 
of attack i H , air-mass flow through the engine Gb, the inlet 
coefficients of flow, resistance, and total pressure (<}>, C x , v g ) , 

Mach number in the inlet opening — Mg, the efficiency of the fuel 
combustion process in the combustor — the combustion efficiency 
in the outlet section and its change along the combustor [illegible] , 
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Fig. 1. 


Diagram of GPVRD. 


the total pressure, temperature, and stagnation coefficients, 
flow Mach number along the combustor duct, the exhaust nozzle 
impulse-loss coefficient, the engine thrust parameters, the 
values of the pressure and friction forces, and other character- 
istics. We shall examine the determination of just some basic 
values. 


We shall set up a characterization of the error in deter- 
mining_parameters by the realtive or absolute standard deviation 
(a or a). Then, if f{xi, X2> . .., is the value under exam- 
ination and depends on the measured parameters X2, • 


<s 






and the value of f falls into the range of values f + Of with 
an error 0.68 relative to o, if the values of a x ^ characterize 
the mean-square random error of the measured parameters . 


Assume we are measuring and Pg' . Then Mjj is determined 
by expression 


_ / rH / i\T» n-* 

' 2 > y-i r l -0(1 + V 7 " 2 * - ' 

and determined under condition 1 ^ ^ " i (which 

is valid at M H >> 1.0), can be found from equation 


«»(< 


□r"£> jT 




I 2 - nl\ 


♦ « 


rl 


(I) 


The incident flow angle of attack ?' H can be calculated using /32 
familiar relations for the degree of static pressure increase in 
the first wedge of the inlet, p = P X /P H - B Y simplifying these 
relations, assuming smallness of weage and shock angles, we can 
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obtain an approximated expression for calculating angle of attack 

v 


l- s 


* - 1 
nT7" 


- O * 


( 2 ) 


where w is the wedge angle of the inlet. Comparison of this 
approximated expression with precise solutions (Fig. 2) revealed 
the possibility of employing it to determine the error of i„. The 
value of cr.£ , following (2) , is calculated using equation: 



Obviously, the error of determining M H has the basic effect on 
a i H * 


The basic inlet parameters <J>, Vg, C x , and Mg are calculated 
using the measured distribution of pressures on the stagnation 
surfaces. For a flat inlet with any number of shocks on the 
stagnation surfaces we can write a system of equations deter- 
mining the flow diagram and the values of <t> and C^. For the 
sake of simplicity and clarity, at = 0 we shall examine the 
case where all shocks other than the first approach under the lip. 
Then 






* 


( 4 ) 


where is the shock slope, oj]_ is the wedge angle of the inlet. 
Using the relation for determining 




n 


and the condition of smallness of angle (i$ «L * «***<. 
we obtain 








( 5 ) 


But since 




then with j 

nTT- 


( 6 ) 


If % 


- x 


I,', it is more correct to determine the absolute error 
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Fig. 2. Comparison of approximated relation (2) 
with precise solution. 

precise solution [6] ; 

equation (2) , this article. 


For analytic determination of the error of the total pressure 
coefficient and Mg, we shall write a system of equations using/34 
gas-dynamic functions (assuming Tg = T£) 


and solving it 


with respect to y^ and A we obtain 

' B g 







Here y(A), f( A), and z (A) are aas-dy^iamic fun 
area of the engine inlet, * « *• ^ 9m 

pressure forces, and R. 


mp 


functions, F BX is the 
ft, are the 
represents the frictioli forces on the 



stagnation surfaces of the air inlet. 


After determining the partial derivatives and simplifying 

the obtained expressions, considering the contribution of each 

term, we obtain the following relations for errors of M and v„: 

g B 




r - 

4 H 

* \[ *♦ JL£- 


(?) 


( 8 ) 


Thus, examination of expressions (5) -(8) provides evidence 
that the basic effect on the error of determining air inlet para- 
meters is exerted by the error in determining M H of the incident 
flow. It is also clear that the errors in pressure and friction 
forces on the stagnation surfaces or have little effect on flow 
parameters in the air inlet opening, since 


1 « 1 . 

Analysis of the errors in parameters determining the effi- 
ciency of operation in the combustor shall be performed using 
the example of combustion efficiency n nc , which is a most impor- 
tant characteristic of operation. 

We should distinguish the peculiarities of methods of 
processing experimental data using a thrust-measurement device 
(permitting measurement or determination by calculation of the 
total flow momentum at the outlet from the combustor 1 ks) and 
without it. 


We shall examine the possibilities of studying DGPVRD in the 
subsonic combustion mode without a thrust-measurement bench. 

The condition Mg 2 = 1 at known air-mass flows Gg and with an 
excess-air coefficient a is necessary and sufficient for calcula- /35 
ting flow parameters in section "G2". In fact, the mass-flow 
equation 



is used to determine temperature Tq 2 and the flow enthalpy J G 
Then, from the energy equation, assuming an amount of heat no 
released upon incomplete combustion of 1 kg of fuel. 


we obtain at a > 1 


*1^ * *«•. (* -*(wc *) at *■ 1 
* Mi*. at • 


T*i * 


(i *®cyc)3 


(9) 

( 10 ) 
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OF POOR QUALITY 


at a < l 




* Hmt. ♦«tL t V3Q l ♦ ^yc)] 


(II) 


Here Q C y C is the amount of heat given off by the engine cycle. 

We shall examine the case a >. 1. Then, from (10) , we obtain 

the formula for the relative error of ri : 

J\C 


*u.*f If- 5 #- *i a J * 

♦ w E I T p . [ 2 , 


“•'<<* «**•* ) ■ 

This formula can be substantially simplified. In fact: 
We can also show that , 


(12) 


-(13) 


Defining 




and considering the results of calculated estimates showing that the 
basic contribution to o^kc i- s ma< ^ e by the error of calculating 
total enthalpy J^/ i.e., the third term in the radicand in (13), 
we obtain an approximated formula for ariKC at a ^ 1: 


®w* 2 {J 
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lI4) 


An analogous relation, but without a in the factor in front of 
the square root, is obtained also for the second case at a < 1. 
We can see from equation (14) that, besides aMjj* the error of 
measuring P also substantially affects the error of ri„_,. 

CjZ J\C 

A more complicated expression for (?n^c is obtained when 
studying the supersonic combustion mode without thrust measure- 
ments. We should stress the need in this case to calculate the 
friction forces along the combustor duct, especially with low 
heat inputs (low engine thrusts) . 
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Questions of determining the error of r) K c in tests with 
thrust measurement, and also questions of comparing both testing 
methods with and without thrust measurement require special experi- 
mental examination and go beyond the framework of this article. 

We shall merely note that use of condition Mq 2 = 1 when processing 
the measurement results in engine tests with thrust measurement 
in subsonic combustion modes introduces an additional error in 
the final result, for example, of B nKc . For instance, we can 
show that in this case the error in calculation of total enthalpy 
would be more than 1.5 times higher than aj£ 2 » determined based 
on a data-processing method not employing condition M<g 2 = 2. 

Examination of the errors in parameters determinig operating 
efficiency in a jet nozzle (for example, the impulse-loss 
coefficient) ini tests of it in an engine system revealed that 
the value of the absolute error is comparable to the value 

of the impulse-loss coefficient itself. This fact demonstrates 
the difficulties of accurate quantitative study of a jet nozzle 
in the engine system. 


Calculations at Mfj = 4-6 were made to determine the level 
of the analyzed errors. It was assumed that the pressure- 
measurement errors in different sections of the flow duct are 
the same and the errors of geometric-value measurement equal zero. 
Friction forces in the combustor and nozzle were ignored in 
calculations of thrust errors in the first approximation. The 
results of the calculations are shown in Figs. 3-6. 



Fig. 3. Relation of M R and o errors to static pressure 
measurement error. 
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The data presented 
illustrate the effect of 
Mh on the errors of 
certain parameters (an in- 
crease of Mh generally 
leads to an increase in 
o ) . The slight reduction 
of the angle-of-attack 
error with increase in Mh /41 
is explained by the more 
flattened path of relation 
(2) at high M H . The effect 
of Mh on the errors of 
air-inlet and combustor 
parameters depends substan- 
tially on the Mh* This 
is because of accumulation 
of errors in flow para- 
meters along the duct. 

For example, if the 
ratio /*« 2 in the 
M H measurement range under 
examination, then at Mh = 6 
in the air inlet the 
analogous, ratios reach 
values 2,9; 

/** - 4.5, a S* /£* =13,5 
(Fig. 4). 

As calculations re- 
vealed, research on a 
combustor in an engine 
system becomes difficult 
when measuring pressures 
with ap = 0.02 + 0.03 
if air-mass flow is not 
measured (Fig. 5) . The 
significant increase of 
the error of g kc with 
_ increase in Mu or a/2 

0 m m 0^ m Sp is explained by reduction 

of the contribution of 

Fig. 4. heat release to total 

energy flow when fuel is 
burned. The cited data 

show that a substantially more accurate value of combustion ef- 
ficiency riKC can ke achieved by measuring the air-mass flow 
through the engine when testing in an attached duct or by prior 
calibrations in engine tests in wind tunnels (errors in other 
parameters are also reduced in this way) . 
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Fig. 5. Dependence of combustion-efficiency error on 
excess-air coefficient at various and a . 

without measuring air-mass flow, 

with measurement of air-mass flow, = 


0 . 02 . 














o afit op* m$ tm 


Fig. 6. Dependence of error in calcu- 
lating internal thrust of a PVRD on 
pressure-measurfement error. 


Possible errors in 
determining engine thrust 
were estimated by inte- 
grating the distribution 
of static pressure along 
the walls of the flow 
duct (Fig. 6) in the 
supersonic combustion 
modes . We mentioned 
above the importance of 
correctly calculating 
friction forces, especi- 
ally at high a or low 
TiKC f i.e., at low values 
of engine thrust. Thus 
the actual or may be 
somewhat higher than 
those shown in Fig. 6. 

Therefore, deter- 
mination of the errors 
in parameters of PVRD 
operating in the 
M = 4-6 range revealed 
that errors in measure- 
ment of pressures in 
the flow and on the 
walls of the flow duct 
have a substantial 
effect on errors in 
operation parameters 
and engine characteris- 
tics. Inaccuracy in 
determining incident- 


flow M H has a definite effect on the values of these errors. 


Our analysis demonstrated that pressure measurement with 
op < 0.02 provides acceptable accuracy of determining PVRD 
parameters and characteristics. 


The authors thank R. I. Kurziner for his useful advice. 
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EFFECT OF IRREGULARITY OF A SUPERSONIC FLOW 
WITH SHOCK WAVES ON FRICTION AND HEAT EXCHANGE 
IN A HYPERSONIC RAMJET ENGINE DUCT 

O. V. Voloshchenko , V. N. Ostras', V. A. Eysmont 


The problem of determining friction and heat exchange in the/43 
duct of a GPVRD [hypersonic ramjet engine] is among the most 
important ones when designing an engine. The flow in the duct 
of such an engine, especially at the beginning of the combustor, 
is substantially irregular along the height and length of the duct. 
The typical calculated pattern of flow in a GPVRD duct in "cold" 
mode is shown in Fig. 1. The pattern was obtained by calculating 
frictionless flow using the method in [1] . As Fig. 1 shows, 
there are shock waves and rarefaction waves all along the duct. 

In the real case this pattern is aggravated by viscous interaction 
and formation of a boundary layer separation zone at the walls 
of the duct. Determination of friction and heat exchange for 
such a flow is an extremely complex task. 

In such situations friction and heat flows are calculated 
with the assumption of flow one-dimensionality, using familiar 
boundary layer calculation methods for supersonic main flows [2] . 

The correctness of this approach when determining friction for 
axially symmetric irregular flows in ducts 7/d ^ 5 long was 
experimentally tested and confirmed in [3]. Bearing in mind 
Reynolds' hydrodynamic analogy, normally used in boundary layer 
calculations [2] , the results in [3] permit the hope that a 
similar procedure is valid for heat exchange as well. However, 
the possibility of employing this determination of friction, and /45 
especially heat exchange, for asymmetric irregular supersonic 
flows, characteristic of GPVRD ducts, required additional testing. 
Furthermore, determination of heat exchange from averaged para- 
meters is known not to provide the distribution of local heat 
flows along the duct walls when shock waves are present, something 
necessary when developing engine cooling systems. 

The purpose of our work was to carry out a calculation study 
of friction and'heat flows in a flat duct with supersonic irregular 
flow, characteristic of GPVRD ducts, considering the detailed 
structure of irregular flow. The possibility of employing fric- 
tion and heat flow calculation from averaged parameters was 
tested simultaneously. 

Selected as the model for investigation was a constant 
cross-section duct with dimensionless length l/Yi =25. A duct 
height h = 0.03 m was assumed as a characteristic dimension. 

The width of the duct was assumed fairly large (b/h = 33) , which 
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made it possible to avoid consideration of the side walls in the 
calculations. Irregularity and asymmetry of flow at the inlet 
to the duct was provided by unilateral deflection of the incident 
flow by a wedge with an angle 8 = 15° . 

The supersonic flow formed in the duct with shock waves and 
rarefaction waves was calculated using the method in [1] . In the 
calculations it was assumed that outside the boundary layer the 
gas is non-heat-conducting, inviscid, and perfect (ae = 1.4), 
there is no interaction between the frictionless flow and the 
boundary layer, and the rise of a boundary layer at the walls 
is compensated for by duct expansion to the thickness of the '.'.boun- 
dary layer displacement. 

The distributions of Mach number M, pressure P, and tempera- 
ture T parameters along the duct walls, obtained through the 
calculations, were taken as the initial data for computing the 
boundary layer. The boundary layer in the duct was considered 
turbulent, formally achieved through appropriate choice of 
pressure at the inlet. The integral method [2] was used to cal- 
culate the turbulent boundary layer. Calculation of the boundary 
layer characteristics with its passage through shocks was per- 
formed using the effective-length method [4] . In distinction to 
[4] , the effective length was determined not from the condition 
of retaining the momentum thickness 5J* = 5J* with the transition 
of the boundary layer across the shock wave, but from the condition 
of retaining friction forces Xfri = Xfr 2 » This procedure, together / 4 6 
with satisfactory agreement with experiment when calculating 
friction in a duct with irregular flow and the more physical idea 
of change in friction force upon transition across the shock waves, 
provides, as demonstrated by calculations, satisfactory agreement 
with experiment [6] as regards the amount of friction stress and 
specific heat flow. The boundary layer was calculated assuming 
T w = const. Calculations were performed for incident flow M ro = 

4+8, to which corresponded a Mach number of an even flow, 
equivalent by momentum and energy expenditure, in the channel 
M = 2.64-4.35. 

Figure 2 shows a typical computation pattern of supersonic 
flow with shock waves and a flow structure in some sections of 
this duct for Moo = 8. As Fig. 2 shows, considerable irregularity 
of pressure P, M number, and velocity direction [illegible] is 
retained throughout this extent of duct. There are four shock 
reflections from each wall in all the modes examined in this 
article with respect to M^, as well as in the mode shown in Fig. 2. 
Estimates revealed that in all the modes the first two shock 
reflections on the upper wall and on the lower produce v 
boundary layer separations. Assessing the effect of separations 
will be discussed below. 
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Weak attenuation of irregularity along the duct is indirectly 
seen in Fig. 3, which shows for M = 8 the distribution at the 
duct walls of relative pressure P/Poo, M, and the parameters 
Pp and Mp of the uniform flow equivalent to it. Figure 3 shows 
that the amplitudes of variations in P and M change fairly weakly 
along the duct. Also shown here are the results of a calculation 
of specific heat flows q w along the duct walls for an irregular 
and uniform flow equivalent to it. 

The results of friction force calculations for each wall 
are presented in Fig. 4 as relations between the_momentum loss 
factor J = [illegible] and duct relative length l = 1 / h and 
Zd = where Xf r = the force of friction, J is the flow 

total momentum in the inlet section of the duct, equal to the 
momentum of a uniform flow equivalent to non-uniform with respect 
to rate, momentum, and energy; d is the hydraulic diameter of 
the duct, for a planar duct equaling d = 2 3h/ [illegible] + h. 

This figure shows the total values of the planar-channel momentum 
loss factors, J = J up + J low . 

We shall first examine the results with respect to friction, /50 
shown in Fig. 4. The distribution of friction forces along the 
duct graphically indicates the effect of shock and rarefaction 
waves on friction force. This effect is strongest at the begin- 
ning of the duct at the points where the shock waves diverge and 
reflect from the duct walls. At these points we find a more in- 
tensive increase in friction force, which then weakens the rare- 
faction wave. A similar pattern occurs in the rest of the duct, 
but with less intensity, a fact associated with weakening and 
damping of shock waves along the duct. A comparison of the 
friction forces at the upper and lower walls reveals that in the 
initial portion of a duct l/h = 2 + 4 long there arises a sub- 
stantial differential of friction forces at the upper and lower 
walls at all the Moo numbers investigated, and this differential 
is then retained all along the duct. The genesis of this differ- 
ential of friction forces at the walls in the inlet section of 
the duct is associated with flow irregularity and asymmetry in 
the inlet section of the duct. Also apparent from the data 
given in Fig. 4 is that the curves of the friction-force distri- 
butions of uniform flows equivalent to non-uniform at 1/ h > 2 
occupy an intermediate place between the values of the friction 
forces for the upper and lower walls. 

There is practical interest in comparing the total momentum 
loss factors in a planar duct with irregular and uniform flows, 
and also with the experimental data in [3] with respect to fric- 
tion in tunnels. This comparison reveals (see Fig. 4) that at 
l/h > 2 the calculated values of friction-related momentum-loss 
factors of an irregular flow differ only slightly from friction- 
related momentum loss of a uniform flow equivalent to irregular 
with respect to rate, momentum, and energy. The obtained calcu- 
lation data also agree satisfactorily with the experimental data 
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in [3] . The slight increase of calculated values over the 
experimental, and also the greater increase in friction along 
the tunnel in [3] , apparently is associated with the presence of 
a boundary layer at the tunnel inlet, increasing in the nozzle, 
and the effect of boundary layer increase in length with 
flow deceleration in the tunnel. 

The above results from calculation of irregular flow 
friction in a planar duct were obtained without considering 
possible boundaty layer separations at the walls. As already 
mentioned, an approximate estimation of the effect of separations 
on the total flow-friction force in the duct was made. Analysis /51 
of the pressure drops in shocks for the irregular flows examined 
in this article revealed that boundary layer separations are 
possible at the upper wall only in the first two shocks, but 
at the lower only in the^first shock. Here the relative length 
of the separation zones Z se p = ^sepl* 1 ' determined from [5], 
may be Z se p - 0.04-0.8. If we assume that there is no friction 
in the separation zones and that the separation zones do not 
affect later flow, but the boundary layer characteristics behind 
the separation zone are calculated similarly to the separation- 
less case, then the reduction of friction through separations 
does not exceed 10%, even with relatively short ducts, as we can 
see from Fig. 4. 

The effect of shock and rarefaction waves on local heat 
flows is graphically apparent when examining the curves of q w 
distribution along the walls of the duct in Fig. 3. Where the 
shock waves fall and reflect from the walls the local specific 
heat flows increase irregularly and exceed by 3-4 times the 
level of heat flows determined from uniform-flow parameters. 

A similar pattern is seen throughout the duct, but with a re- 
duction in the maximum level of local heat flow, associated with 
damping of the shock wave intensities. The difference in 
distribution of specific heat flows at the upper and lower 
wall of the duct is associated with irregularity and asymmetry 
of flow in the inlet section. As a whole, the distribution of 
specific heat flows at the duct walls reflects the nature of 
pressure distribution by length in an irregular flow and can be 
approximated by relation 


kt-i * ^ Pi !>c-i ^ * 

where (P^/P^_i) is the local pressure increase in the shock 
at the duct wall. 

Figure 5 shows a comparison of the calculated change of 
heat flows in the shock waves of this duct with the experimental 
correlation for separation flow q-^/q^-i = (P^/P-i-i) n . . . [6], where 

the exponent n for the point of separation "attachment" equals 
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0.85, and for the "plateau" region, n = 0.625. Also shown are 
the calculated points obtained for the conditions of an 
experiment [6] in which heat flows were measured with = 8.6 

of a flow over an internal obtuse angle. 

As Fig. 5. shows, calculation of local heat flows for a 
frictionless flow produces intermediate values for heat flows, 
obtained with separation (between the "plateau" and the /54 

"attachment" point) . The obtained level of calculated local 
heat flows is 20-40% below the maximum heat flows at the separa- 
tion attachment point. 

Figure 6 shows, for = 6 and 8, the ratio of the 
integral heat flows for each duct wall along its length for an 
irregular flow and the uniform flow equivalent to it, Q = 

(/<3w^ J/^qwpd • We can see that in the initial section of the 
duct (l/h = 2-7) there is a substantial increase of integral 
heat flow at the upper wall over that of the lower and uniform 
flow, and it diminishes at length l/h ^ 8. Here the integral 
heat flows at the lower wall throughout the duct remain some- 
what below those at the lower and those calculated using one- 
dimensional parameters. The total heat flows along the duct, 
calculated for irregular flow and the uniform flow equivalent to 
it with respect to rate, momentum, and energy, coincide approxi- 
mately at length l/h .> 4. 

Thus, as research demonstrated, at l/h ^ 4 calculation of 
total heat flows and friction forces in the duct of a GPVRD 
can be made using the assumption of flow one-dimensionality in 
the duct and familiar boundary-layer computation methods. This 
approach produces satisfactory results for both axially symmetric 
and planar asymmetric flows. 

The local specific heat flows should be determined using 
local flow parameters, since their level may be 3-4 times higher 
than that of specific heat flows calculated using uniform-flow 
parameters . 
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FEATURES IN THE SELECTION OF BASIC PARAMETERS OF 
GAS-TURBINE-ENGINE COOLED TURBINES 

O. N. Yemin 


1 . 


Peculiarities in calculation of cooled turbines 


1 


/56 


The extensive use of evolved air-cooling of the flow 
portion of the gas turbines of aviation GTD [gas-turbine engines] 
has necessitated a more precise procedure for selecting their 
optimal parameters. As applied to uncooled turbines, parameter 
optimizations have been fairly thoroughly developed and are 
discussed in scientific-technical and academic literature [2, 3]. 
Works [4-6] have cited data needed for preliminary calculation 
of a cooled turbine or a turbine with cooled stages. 


In addition, the initial results from use of cooled turbines 
indicated a need for considering in the choice of their basic 
parameters the features of their operation caused by the presence 
of a system to cool the flow-area components. 


Calculations demonstrated, for instance, the need for more 
accurate assessment of the cooled-turbine r.p.m.'s adopted 
for the design. 


Therefore, this article examines research at the engine 
department of MAI [Moscow Aviation Institute] concerning the /57 
effect of r.p.m.’s on the efficiency of a cooled one-stage turbine. 

The analysis was performed for a one-dimensional and the 
simplest quasi-two-dimensional model of the turbine's operation, 
i.e., considering its parameters at the mean diameter. The 
feature of this outline of calculating a cooled turbine using 
the mean diameter consists in initial calculation of an uncooled 
stage with determination of the velocity-triangle parameters 
at the mean diameter, the diametric and fractional (axial) 
dimensions of the flow-through portion (height and width of the 
blddes, widening angles of the flow-through portion, etc.), and 
then, for the thusly calculated "uncooled analog" of the turbine 
being designed, determination of the cooled air flow required to 
obtain the turbine temperature with which the given long-term 
strength margins are achieved. Later the working-medium flows 


^"The development of aviation and cosmonautics is characterized by 
a constant increase in parameters of powerplant cycles, especially 
temperatures in the main combustion chambers of gas-turbine 
engines [1] . 
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are ascertained in the calculated sections of the stage, in 
accordance with the cooled-air delivery outline, blade height, 
and then velocity-triangle parameters. 


Thus, although calculation of a cooled turbine according to 
mean diameter is reduced basically to calculating an uncooled 
turbine with additions and adjustments, the values of the optimal 
(expedient) parameters selected when designing a cooled turbine 
may be substantially different from the parameters of the 
"analogous" uncooled stage ("uncooled analog"). By "analogous" 
uncooled stage we shall mean a stage with the same values of 
basic structural and gas-dynamic parameters (H<i>u/ PT' a l» $2' 
u, Dgp, h n*..), but made uncooled (for example, of a hypo- 
thetical superheat-resistant material) . 

Such a representation is extremely arbitrary, since a cooled 
turbine, for example, despite the presence of a cooling system 
(equipment for delivering and distributing cooled air) , should 
be viewed as geometrically similar to an "uncooled" one in which, 
of course, this equipment is not provided. As will be shown 
below, however, this comparison is useful, since the efficiencies 
of uncooled stages and turbines are determined fairly firmly 
and are well-confirmed by many years of experience in their 
construction and use. The relation shown in Fig. 1 is used in 
our analysis [4] . 

Determination of a cooled turbine's efficiency is difficult, 
especially in the initial design stage, since it naturally depends 
on the peculiarities of the cooling system, which are yet unknown. /5 9 

Thus, in the first approximation it is suggested that a 
cooled turbine's efficiency be defined as 


n; 


= n* 


(l) 


'T.c T.unc c 
where 0*T.unc i s the efficiency of an uncooled turbine stage 
with basic parameters (Hi> U / Pt> the same as the cooled 

one under consideration, with the same nozzle-arrangement and 
impeller geometry, called the "analogous uncooled" stage above; 
ri c is the relative efficiency, considering all the peculiarities 
of the cooled turbine's operation (inequality of flows in 
stage components, temperature reduction with gas and cooled-air 
mixing, additional losses with mixing, etc.). 


To determine ri c , it is proposed that the experimental 
(statistical) relation shown in Fig. 2 be employed. Obtained 
on the basis of processing a relatively small amount of experi- 
mental material for the most typical outline of convective-film 
cooling, it is approximate, but adequate for analyzing turbine 
parameters at the stage of engine draft design. The value of 
ri c depends chiefly on the relative total flow of cooled air — 

G y (to the stage) - — so the value of n is sometimes arbitrarily 

C L C 
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Fig. 1. Efficiency of an "uncooled" turbine 
stage. The broken line shows the graphic 
determination of the ratio between r.p.m. 
and efficiency at maximum setting (A) and 
in the maximum efficiency mode (B) . 


represented as the ratio of a turbine with the given total 
efficiency value to that of the same turbine with G c j~ = 0. 



G c.CA + C c.PK 
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The required cooled-air flow can be found in the very first 
step of determining the basic turbine parameters using the familiar 
relation, the equation associating the coefficient of temperature 
efficiency of cooling 0 with the relative air flow needed to cool 


one blade [5] : 





Fig. 2. Effect of total relative amount of 
cooled air G c .cA + G c> px at a stage with convec- 
tive film cooling of blading on its relative 
efficiency . 


2, Determination of Optimal (Expedient) Turbine Parameters 
Providing High Efficiencies 


Optimization of parameters in accordance with the employed 
calculation outline assumes obtaining the greatest value of 
parameter [symbol omitted], consisting of two factors, each of /61 
which depends on thrbine parameters associated by extremely 
complicated relations. 


The efficiency of an uncooled turbine depends chiefly on 
the dimensionless parameters describing the process in a simple 
stage at the mean diameter relative to the height of the flow- 
through portion and the law of flow-portion profiling by radius. 




’• unc 


s i ^ TU * ft » 




(3) 


For the given value of theoretical work Htu of gas flow G 
and the given or assumed value of the turbine's mean diameter, 
the choice of some definite value of the flow coefficient C a 
(or its associated value of angle ot]_) determines the value of 
Dcp/h£, and thus also the degree of reactivity at the calculated 
mean diameter. In this case the efficiency would be determined 
by the loading of the stage, i.e., the value of the peripheral 
velocity or, consequently, the turbine rotor r.p.m. 
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Reduction of stage loading (i.e., increasing the r.p.m. 
at a given mean diameter value) usually leads to a decrease in 
the efficiency of an "uncooled" turbine stage (see Fig. 1) . 

The tensile stress at the root of an operating blade increases 
because of centrifugal forces [2] . 

e, « ^ (4)’ 

2 

where is the annular area at the stage outlet, m ; $ is the 
form factor of the operating blades ($ = 0.58-0.62). 

Obviously, to provide the required strength margin for a 
turbine with a higher r.p.m. we have to ensure fairly high values 
for the long-term strength of blade material. Of course, high 
long-term strength values can be provided by reducing the tempera- 
ture of the blade material, and this requires increasing the 

cooling air flow. This causes a reduction of fT . 

c 

Recall that the value of n* T>unc increases when r.p.m. 
increase. 

Thus, to obtain a high efficiency value for a cooled turbine 
the r.p.m. should be chosen so as to produce the highest product 
ri^.unc’^c' consisting of values changing in opposite ways with /62 
change in n>r. 

Generally speaking, solution of the problem remains diffi- 
cult. Thus, in each specific case it is expedient to perform a 
series of variant calculations. 

A typical example of selecting the optimal r.p.m. value is 
analyzed in detail below. 

3. The Concept of Limiting R.P.M. 


As an example of the procedure described in [6] , calculations 
were performed for a gas-generator turbine intended for creating 
a TRDD [ducted-fan turbojet engine] with a total pressure ratio 
in start conditions = 30 and a combustor temperature 

TjS = 1700K. The pressure ratio in the gas-generator compressor 
^KrT = 11*5 and the gas flow Gp = 65 kg/sec. The convective- 
film cooled blades are made of alloy ZhS-6, the cooled-air tempera- 
ture is assumed T£ = T]£ = 880K, and the strength margin K<jp = 

2 for a calculated operating time r = 250 hrs . The mean diameter 

of the turbine is assumed D = 0.7 m. 

cp 

We note first of all that experience in performing systema- 
tized variant calculations reveals the expediency of varying 
both the value of the turbine mean diameter and the value of the 
blade-material temperature, T £. The temperature of the diaphragm 
vanes was assumed T^ = 1300K in all calculations.) 
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Fig. 3. For defining the idea of limiting r.p.m. 

of a turbine — = const. 

a) input velocity triangle element; 

b) calculation of operating impeller-blade strength 
margins, under the condition that [expression 
omitted] . 

As Fig. 3 shows, to meet condition K(jp = 2 at the assumed /62 
blade material temperature and with the given gas flow value 
through the stage, increasing the r.p.m. makes it necessary to 
increase angle a, and thus to decrease blade height. When 
angle aq is increased, however, the velocity in the impeller 
clearance increases, as does its axial component. When the 
limiting expanding power of the nozzle equipment is reached 
(axial component equal to the speed of sound) , further reduction 
of the flow-section area, as shown in Fig. 3a, and thus also 
of the blade height, is curtailed. Thus, for highly loaded 
turbines, in which the value of Cq u is significant, and especially 
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for gas-generator turbines, the appearance of limiting r.p.m. 
becomes characteristic, since there is no need to compensate 
for further increase in r.p.m. by reducing blade height. In 
other words, limited nozzle-capacity modes may become charac- 
teristic of these turbines, and with these it is obvious that 
the turbine diameters have to be increased when r.p.m. are 
increased. The limiting value of angle aj_, determined by the 
formula following directly from Fig. 3a 


tan* 


♦lim 



(7) 


Corresponding to the limiting r.p.m. value are the effi- 
ciency value points limiting for:. the given temperature. By 
connecting the limiting-efficiency points for various temperatures 
we obtain a curve of limiting efficiencies (Fig. 4) . The relation 
between n^lim and has a distinct maximum. The presence of a 
maximum is explained by the opposite effect on efficiency of 
the loading factor and the amount of cooled air. The nature of 
the dependence of these values on r.p.m. in limiting modes is 
shown at the top in Fig. 4. At the same time, for given values 
of T i the dependences of limiting efficiency on r.p.m. may 
have an extremum, and we have to determine, in addition to the 
limiting efficiencies, the values of the maximum efficiencies, 
which generally may differ from the limiting. 

4. Relation Between Maximum and Limiting Efficiencies 


As we already mentioned, at T i - const, an increase in 
angle (to as much as oti.iim) leads to a monotonic increase in 
r.p.m., corresponding to the assumed strength margin K = 2.0, 
right up to the value n = n^^. y 

As we know, however, the turbine-stage efficiency has an 
extreme character when the calculated value of angle a is changed 
[2] . The position of this extremum is determined by the config- 
uration of the velocity triangle, which determines the amount 
of loss in the stage and the exit velocity (C2 and X 2 ) . The 
location of the extremum of function = /(«]_) is not associated 
with the strength indices of the projected stage and the value 
of ctl opt (by efficiency) can be found "to the left" and 
"to the right" of the value of »i.iim (by strength), as shown 
in Fig. 5. A relatively simple method of establishing the 
relation between maximum and limiting efficiencies is provided by 
the simplified procedure of evaluating the efficiency of a cooled 
turbine, in accordance with which the efficiency of an "uncooled" 
turbine n$. unc is determined using the graphic relation shown in /67 
Fig. 2. 

Preliminarily, we establish that along line Tjr = const, in 
Fig. 4 we can consider valid the following approximated relation: 
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Fig. 4. Explanation of extreme character of 
relation between cooled-turbine efficiency and 
r.p.m. (a) . 'Possible relationship of locations 
of its maximum and limiting efficiencies (b) . 


Numbers indicate lines of constant values of 

temperature of blade material (T., > T-,„ >... T 7r ). 

b b 2. bo 


* L ; Gc 




= > ct>ni , 


The assumption that the relative flow of cooled air remains /67 

unchanged is, generally speaking, approximate, since change in 

r.p.m., and thus in turbine circumferential velocity u, 

strictly speaking, leads to a change in the stagnation temperature 

in relative movement — T* . 

w 

In addition, for the assumed condition K 0 p = 2.0 the value 
of the tensile stresses at the blade root should also be un- 
changed: 

y ^ x . 


Converting the simplest formula for determining stresses at 
the root of an operating blade, we obtain 







* h 

" 5 £ warn* 
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Fig. 5. 
case 1 , 
case 2, 


Relation between n* 


and 


W ^ ere nr i® ax n lim an ^ h'fmax 
where n ^ < 

n*max n, and 


Tmax 


> 


n Tlim 


lim 


Thence it fol- 
lows that at T i = 
const, we observe 
condition 

C. . (8) 

On the other 
hand, when designing 
a turbine with a 
certain value of 
theoretical work 
Ht u and a given 
mean diameter D C p, 
when the r.p.m. 
is varied we observe 
condition 




( 9 ) 


Combined solu- 
* tion of equations (9) 

n Tlim and (8) leads to 

= ^Tlim equation: 


*** K ' (id) 


By applying this relation to the graph of the efficiency of an 
"uncooled" turbine (see Fig. 1) we obtain the possibility of 
establishing the relation between r.p.m. in limiting mode (point 
A) and the maximum efficiency mode (point B) . 

On this graph, point B defines the point where the curve 
H Tu Ca = const. touches one line of constant values of n^.nnc' 
and at the point of tangency (point B) the maximum value of 
ri£ unc is achieved. Obviously, the maximum value of the effi- 
ciency of a cooled turbine r)$ >c would also correspond to this 
mode. This follows from the fact that all the examination is /68 
made along line T i = const., for which, as shown above, the 
relative cooled-air flow G c and the relative efficiency, consi- 
dering the effect of cooling, fi c , remain unchanged. 

Thus, with the same values of T £, between the parameters 
of the cooled and uncooled turbines there is an association 
defined by relations 
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(II) ; 

To determine a rational value of the turbine rotor r.p.m. 
of a gas generator, we have to construct the relations of 
maximum and limiting efficiencies to r.p.m. Using relations (11) , 
we note that for some combination of turbine strength indices 
and relations between its efficiency and r,p.m. the line 
T i = const, may be in the form of the broken line in Fig. 4, 
i.e., with an extremum at point C, lying below pount A on the 
limiting efficiency curve. Calculations done according to this 
procedure for a gas generator with the parameters given in 
section 3 revealed that, to obtain the highest efficiency value, 
we would have to select the following values of the operating- 
blade material temperature, turbine r.p.m., stage loading factor, 
and total relative cooled-air flow: 

1200K 

13.000 r.p.m. 

2.0 

4.5 

This would provide a turbine cooled-stage efficiency of 

0* = 0.865. 

T.c 


T = 
n = 

Htu = 

G = 
c 


It 

r t T 1 im 
n 

n. . 

lim 


It unc ft 

1 * unc * 

Ca. 


However, this question of choosing expedient parameters 
for a gas-generator turbine cannot be considered answered, 
since we have to introduce into the examination, in addition to 
turbine efficiency, the compressor operation indices, also deter- 
mined by the chosen r.p.m. 

We can also assume that the expedient r.p.m. would be found 
near the extreme values of the curves of TVf.iim and hip^x/ which /69 
are generally quite close to each other. 

The above sample calculation allows formulating some con- 
cepts about the features of the basic parameters of a high- 
temperature gas-generator turbine stage. 

The calculations demonstrated convincingly that the internal 
efficiency of a cooled turbine, as opposed to an uncooled turbine, 
is affected basically by, in addition to the loading factor, 
the amount of cooled air blowing on the surface of the blade, 
so the attempt should be made to increase the blade's design 
temperature and reduce the percentage of discharged air. For 
reducing G c we can sometimes recommend measures to decrease the 
cooling-air temperature in a special heat-exchanger [ 6 J . 
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The highest cooled-turbine efficiency is achieved in modes 
with which the value of the loading factor is higher than that 
recommended for uncooled turbines . 

The minimum value of turbine-blade height for the given 

temperature of material corresponds to the limiting curve. But 

in maximum efficiency modes, too, high-temperature turbines of 

gas generators are characterized by a high value of ratio 

D /h 7 = 15-25. 

cp l 

The extreme nature of the relationship between efficiency 
and r.p.m. leads to the need for designating a "tie" of the engine 
rotor to determination of the expedient rotor r.p.m. of its gas- 
generator . 
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ON SELECTING THE OPTIMUM TOTAL WEDGE ANGLE 
FOR A HYPERSONIC RAMJET ENGINE AIR INLET 

V. I. Penzin 


The most important component of a ramjet engine (PVRD) , /70 

determining its efficiency, is the air inlet. At high flight 
speeds, with which it becomes expedient to employ supersonic 
combustion in the combustor (PVRD) , inlet requirements change 
markedly, since the optimal inlet shape becomes dependent on the 
engine's operating conditions, i.e., on the excess-fuel ratio 
(stoichiometric equivalent $) [1]. A change in the air inlet's 

shape entails both a change in the combustor shape and in the 
level of total pressure loss in the combustor [2] . In [1] an 
examination was made of a GPVRD [hypersonic ramjet engine] with 
the simplest two-dimensional, one- and two-shock air inlet. 

This article investigates the effect of stoichiometric equi- 
valent $, number of shocks N, relative inlet throat area F r , 
lip shape, and other parameters in the Moo = 6-10 range on the 
optimum total wedge angle of a multishock air inlet 0 R . 

Figure 1 shows a computational outline of a GPVRD with a 
two-dimensional multishock fixed air inlet. Taken as the optimum 
value of 0 k was the value at which the thrust coefficient C a 
is maximal at the given $ . The GPVRD thrust calculations 
considered friction pressure forces on the stagnation surface 
and internal and external surfaces of the lip. All of these 
forces are easily calculated for two-dimensional flow. The 
internal and effective (less external resistance) engine thrusts /72 
were determined with the assumption of one-dimensional flow in 
the combustor and nozzle. 


The initial throat area (combustor) of the air inlet was 
determined without considering the effect of viscosity, assuming 
the speed of sound in the throat upon start. Then the value of 
^r.init was increased to calculate the effect of viscosity [3] 
and to provide inlet operation in the M range, after which the 


effect of this change on 


K.opt 


was determined. 


Hydrogen was examined as the fuel, the available heat factor 
(complete combustion) was assumed equal to 0.9, and the nozzle 
momentum factor n n = 0.95. The optimal shape of the GPVRD 
flow area (air inlet and combustor) does not depend on the area 
of the nozzle outlet section, and in the calculations the latter 
was assumed equal to the area of the entrance to the inlet. 


An increase in inlet angle 0^ leads to increase in the 
intensity of shock waves and loss of total pressure in them, but 
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heat losses in the combustor decrease, since the flow velocity 
in the chamber before the fuel-delivery section diminishes. 

At 0£ = 0 k. opt the total pressure losses in shocks, and with 
absorption, reach a minimum, and thrust reaches a maximum. 

Figure 1 shows the relation between 0 k. opt and stoichio- 
metric equivalent $ at M m = 10, using an air inlet with 1 to 
infinity shocks. In the latter case we examined an" isoentropic 
air inlet (inverse of Prandtl-Meyer flow) . The graphs shown in 
Fig. 1 and Fig. 2 correspond to the condition of equal intensity 
of shocks and are calculated for two cases: 1) the difference 
between the total wedge angle 0 r and the undercut of the lip 
are constant (0^ - 6 = 0.16 radians) (solid lines); and 2) 
zero undercut of lip, 6=0 (broken lines) . An examination of 
Fig. 1 reveals that the number of shocks in a multishock air 
intake substantially alters the optimal value of 0^. For example, 
a move from N = 1 to N = 4 changes 0 k. opt at $ = 1 by A0 K - 0.5 
radians. Reduction of the lip undercut leads to an increase in 
the value of ©K.opt* But it should be stressed that condition 
6=0 (at N > 2) is apparently impossible to bring about in 
practice because of the high value of 0 K - 6 , which may lead to 
stall. Increasing the heat input (increasing $) leads to an /74 

increase in 9 k. opt' and more so as 6 is smaller. When $ is changed 
from 0.5 to 1.6 the value of 9 k. opt ma Y increase (at N = 3) 
by 0.1-0.15 radians. 

At = 10 the value of 0 k. opt is reached at practically all 
values of $ <. 1 . At M ro = 6 (Fig. 2) and with 0^ = 9 k. opt i n a 
constant-section duct, only a relatively small amount or fuel 
can be delivered ( Him < 0.2-0. 3) . Duct choking occurs when $ 
is increased above the limiting value, and the area of the inlet 
throat and the flow velocity in it increase. The left branches 
of the curves in Fig. 2 correspond to the optimal value of 0 K , 
and the right to non-optimal values (choking mode) . The conclu- 
sions drawn from the results of calculations at M =10 can be 

oo 

extended to the case where = 6 if we examine only the left 
branches of the obtained curves. Additional information in Fig. 

2 relates to the effect of lip wedge angle y (see Fig. 1) . In 
the calculations corresponding to Fig. 2 it was assumed that the 
length of the lip segment with undercut l Q j-, equals 0.15 times 
the height of the entrance to the inlet. 

It is apparent from Fig. 2 that change in y from 0.07 to 
0.27 radians leads at $ = 0 . 2 to an approximately 0.15-radian 
reduction in ©K.opt* it is also interesting to note that the 
right branches of relations 0 k = 0k ($) at N > 3 are close to 
each other and can be replaced by one curve for preliminary 
evaluations . 

Generally speaking, it is inexpedient to increase $ > $ii m 
by reducing 0 k in a GPVRD, since this reduces engine efficiency. 
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It is more proper to accomplish heat input in an expanding com- 
bustor after section F = const. In this case, at $ > the 

value of 0v would be constant and equal to (6 k. opt) max* In 
practice, however, cases are possible where increased heat input 
into an F = const, chamber may be useful. This may be because 
of the advantages of simplifying the fuel-delivery system, 
possible increase in combustion efficiency, etc. 

Figure 3 shows generalized relations between 9 k. opt and 
for cases N = 2 and 3 (broken and solid lines) for various values 
of stoichiometric equivalent $. The envelope of the maximum 
values of ^K.opt dirked off by the broken/dotted line. An 
examination of the graph makes it apparent that 0K.max is reached 
at = 9, and it diminishes with further in crease in M , a fact 
explained by the reduction of relative heat input (under the con- 
dition that $ £ 1) . 

Optimization of a GPVRD-powered cruise aircraft at = 8, /76 

10, and 12 is conducted in [4], A two-shock, two-dimensional 
air inlet with zero lip undercut was examined. The optimal 
values of stoichiometric equivalent $ and the total lip wedge 
angle 9 r at 8 and 10 were obtained as equalling 0.45 and 0.53, 
0.295 and 0.23 radians, respectively. These values are marked 
out by circles in Fig. 3, and thence we can see that the values 
of 0K.opt from [4] are somewhat lower (by about 0.05) than those 
in our work. This can be partially explained by the fact that 
in [4] the air inlet was assumed variable at start, i.e., with 
a smaller throat area. 

In actual air-inlet constructions, the condition of shock- 
wave equal intensity is not strictly observed. The slope of the 
first panel of the deceleration wedge is made larger than 

is dictated by the condition of equal shock intensity, and the 

air inlet becomes shorter. Figure 4 shows relations between 
C a and the slope of the first panel of a three-shock air inlet 
0^. The slopes of the second and third panels were selected 
from the condition of equal intensity of the second and third 
shock waves. A strong dependence of C a on 0p can be seen, 
when 0^ is smaller than the optimal value. In the limiting case, 

at 0p = 0, the air inlet degenerates into a two-shock one. At 

®1 > ®1 opt t ^ ie va lue of C a changes slightly. Also in Fig. 4, 
the broken/ dotted line marks out the relation between 0 k. opt and 
0]_ at $ = 1. We can see from the graph that an increase in 0^ 
leads to an increase in 9k. opt- A change in 0]_ from 0.1 to 0.2 
radians leads to a 0.07 increase in 0 k. opt- The preceding 
calculations were performed in the assumption that the inlet 
throat area is determined under the condition that M in the 
throat equals unity upon start. The throat area would be greater 
than initial also in the case where an invariable inlet is used 
in the M range. 

CO 
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Figure 5 shows the relations between 9 k. opt a ^d relative 
throat area of three- and four-shock air inlets at various values 
of $. An increase in the throat area leads to a slight increase 
in opt* since the actual difference between F r and the initial 
value is significantly smaller than that shown in Fig. 5. For 
example, the actual air-inlet throat area in the condition of 
turbulent flow in the boundary layer exceeds the initial value by 
approximately 8% [3] . The need for air inlet operation at 
Mco = M [symbol omitted] -3 leads to an approximately 20-30% 
increase in F r . Considering that the slope of the curves of re- /80 
lation 0 K = e K(Fr/F r . init) is practically independent of $, 
a 30% increase in F r leads to a 0.01-radian increase in 9 k. opt* 
Because of the small change in the thrust coefficient C a near 
the maximum with change in ©k, when determining 9 k. opt i n para- 
metric calculations we can use the initial value of the air 
inlet throat. 

In conclusion, we note the effect of the shape of the fore 
part of the lip on ©K.opt* Figure 6 shows relations between 
©K.opt an< 3 ^ f° r various values of T 0 j-, and y at M =10 and 
N = 3. Change from the shortest and ^thinnest" (Job = 9.05 and 
y = 0.07 radians) to the longest and "thickest" (Z Q b = 0.15 and 
y = 0.14) fore part of the lip at $ = 1 leads to a 0.07-radian 

reduction of 9 . 

J\ 
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THE IDEAS OF F. A. TSANDER AND ASSESSING THE USE 
OF JET ENGINES FOR ACCELERATING AEROSPACE AIRCRAFT 

Yu. M. Annushkin 


1. Basic Ideas of F. A. Tsander. and the Future of Aerospace 
Science - 1 - " 

A basic idea that can be ascribed to F. A. Tsander is the /81 
suggestion that atmospheric air be used in jet engines to accel- 
erate aerospace aircraft. His design [1] published in the 1920 's 
is widely known today. 

Of course, the level of development of aviation and rocket 
technology at that time made it impossible to carry out space 
flights with a rocket or an aircraft launch. The creation of 
space-rocket systems became possible only in the 1950's, and 
creation of more economical winged aerospace booster systems is 
apparently a matter of the near future [2] . 

In the 1930's reaction-engine technology was split into 
two distinct fields: the first (use of VRD [jet engines] and 
winged craft) was assigned to aviation, the second (use of ZhRD 
[liquid-fuel rocket engines]) merged into independent rocket, 
and later aerospace technology. The problems of creating rockets 
with ZhRD in a system were simpler than those of creating 
vehicles using the Earth's atmosphere to produce VKS [aerospace 
plane] lift and VRD thrust. A significant role in this was played /82 
by the fact that ballistic missiles were an even more important 
weapon than winged craft. 

The development of aviation and aerospace science and tech- 
nology is associated most of all with the development of aviation 
and rocket engine-building. Understanding this, Tsander, especi- 
ally in the last years of his life, focused his attention on de- 
velopments of 'and research on jet engines which could be employed 
to boost spacecraft. The status of aircraft engine-building in 
the 1920 's did not allow him in his designs to examine engines 
using atmospheric air at speeds exceeding 400 m/sec. This appar- 
ently explains the fact that in his design for an air-launched 
vehicle he adopted propeller-driving powerplants [VMG] employing 
air. 


Materials from Yu. S. Voronkov, T. M. Mel'kumov, Yu. V. Biryu- 
kov, and I. A. Merkulov, published in collection [3], were used 
in the preparation of this section. 
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Aircraft engine-building theory and practice have now 
entered a region where we can examine the use of VRD, especially 
hydrogen VS [boost stages] , in aerospace planes [VKS] 
and aerospace transport systems [TAKSi] at flight Mach numbers 
Mf = 10-12 [4, 5]. Present development of space-rocket technology 
is taking the road of designing cheaper recoverable spacecraft 
[KAMP] , in the first generation traditionally employing ZhRD, 
in the second, VRD (TRD [turbojet engines] , PVRD [ramjet engines] , 
and GPVRD [hypersonic ramjet engines] ) . 

Aerospace planes with hydrogen VRD as part of their SU 
[powerplants] may be single-stage with an air launch and airfield 
landing after returning from space flight; TAKSi may accelerate 
a rocket to velocities V = 3000-3600 m/sec and return "to base." 

Use of aircraft designs with VRD in space-rocket technology 
is, at the current stage, dictated considerably by the need con- 
duct space research for assorted economic tasks. In this connec- 
tion it is becoming important to carry out space flights that 
provide proper aircraft "comfort" for the crews and expand the 
range of "launch points." It is for this reason that work is 
under way to create winged booster rockets and aerospace planes, 
i.e., self-launched and -landing aircraft. These operations are, 
in principle, a continuation of the ideas of Tsander, who first 
pointed out the substantial energy and operational advantages of 
using an aircraft over a rocket launch, suggested and thoroughly 
substantiated by K. E. Tsiolkovskiy . It is even more significant /83 
that Tsander proposed using Aerodynamic braking of an air- 
transport craft when it returns from space flight (without starting 
the rocket engine) . Tsiolkovskiy and Kondratyuk later re j 
ferred to this idea many times [3, p. 28]. 

Hydrogen-powered VRD (or combined RD {.jet engines] with 
hydrocarbon and boron hydride fuel for VKS) using ZhRD in the last 
stages are a new advance in the development of aerospace tech- 
nology [5] . 

Without making an absolute out of experimentally obtained 
knowledge, we should note that it plays a definitive role in 
technology, especially aerospace technology, so practical 
creation of VKS and TAKSi will be possible only after prolonged 
and difficult basic and experimental research on test-beds and 
in flight, particularly on model VRD and VKS. Questions of 
providing the aerospace industry with scientific and practical 
projects are being given paramount importance throughout the 
world. Among the most important problems in engine design is 
that of setting up highly efficient fuel combustion in engines, 
with minimal longitudinal dimensions and provision of heat- 
resistance; also important is the porblem of cooling the airframe 
structure with fuel (hydrogen) during flights at high supersonic 
and hypersonic velocities (M^ = 3-12) within the atmosphere. 
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The problem of organizing the process in hydrogen PVRD 
combustors and providing their high economy is examined below. 

2. Hydrogen Ramjet Engines for Boosting VKS and TAKSi 

a) Peculiarities of VRD for VKS 

Calculations for evaluating the economy of different types 
Qf engines relative to flight speed, taken from [5] , reveal that 
VRD are much more economical than ZhRD (Fig. la) . For example, 
use of a combination of a kerosene-powered TRD with a hydrogen 
PVRD, shifting to GPVRD mode at Mf - 6, allows obtaining a ZhRD 
specific pulse 10-4 times higher than that of a ZhRD in the 
M = 0-12 range, respectively. Use of only high boiling-point 
components (of the kerosene or boron hydride type) as fuel for 
VRD also substantially increases the economy of the engine over 
a ZhRD, but it should be remembered that at Mf > 4 these fuels, /85 
because of their low cold-margin, are incapable of ensuring 
operation of VKS engine structures and components to high Mf . 

Up to Mf ^ 5, however, use of such fuels is possible, and abroad 
they are viewed, with sufficient foundation, as the probable 
fuels for winged aircraft having various purposes. This is all 
the more attractive a proposition in that such use eliminates 
the problem of storing winged aircraft and SU in the fueled 
condition . 



* 2 s * ® n„ 


Fig. la. 
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Aerospace planes are large in volume and area, in terms of /85 
their good aerodynamics, especially for high-altitude flights. 

They can be refueled with liquid hydrogen immediately before 
launch and there is no need to keep VKS in long-term storage 
in the fueled condition. Hydrogen is not the only fuel capable 
of assuring powerplant and airframe component cooling in hyper- 
sonic flight inside the atmosphere. 

Combined and ramjet VRD for VKS and TAKSi may be extremely 
large and integrated with the aircraft's surfaces for their maxi- 
mum use to compress air ahead of the intake and expand the gas 
behind the engine nozzle (see Fig. lb). 

To provide possibilities of experimental development on 
test-beds, the powerplant may consist of identical linear models, 
one of which (or a reduced-size, geometrically similar model) 
can undergo test-bed development or be studied on a flying 
laboratory (as is proposed for conducting research on two- 
dimensional GPVRD models in the U.S.A. on aircraft of the 
X-15 and X-24 type, see [5]). 

The modular construction of VRD is very convenient for test- 
bed research, and as experience has shown, it allows obtaining 
reliable scientific data at low cost. 

Discovery of the main criteria determining the efficiency 
of the process in VRD is at the top of the agenda for scientific 
experimentation; it is vital that modeling criteria be properly 
chosen and that experience in model research be employed in 
designing full-scale VRD. It is more important to discover 
the modeling principles for VRD than for other types of engines, 
since their thrust characteristics depend on a large number of 
factors [9] , especially flight conditions (Pg, T^, Mf ) , air- 
intake characteristics (<f>, max) » type of fuel (H u , Lq) , /87 

relationships between air ana fuel consumption (a = a B /a T LQ 
or 3 = 1/a), relative flow-duct areas (see Fig. 2) (F C /F H , 

Fjq/Fr, F cr /F B , F cr /F c , etc.), fuel combustion efficiency in the 
combustor, ^r e f' which is the main characteristic of the opera- 
ting process. 

The relative combustion efficiency in the engine (nr e f) is 
determined by the thrust, measured in test-bed tests, and consi- 
ders, in addition to chemical underburning, losses associated with 
irregular parameters at the nozzle exit and heat losses in the 
combustor, and it generally depends on almost all the afore- 
mentioned factors, especially a and combustor construction 
(number of nozzles in the fuel-nozzle unit, uniformity of their 
distribution by chamber section, chamber length, area, etc.). 

This specific feature is usually not accounted for in general 
VRD theory. Existing methods of calculating VRD parameters and 
thrust characteristics are either based on the presumption of 
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plane 


Fig. lb. 







an ideal process in PVRD combustors (nfp,f = 1) , and provide for /87 
determining maximum possible thrust characteristics, or with 
parametric variation of ripef make it possible to discover the 
effect of nref on VRD thrust characteristics (see [6] , [7] , [8] , 

for instance) . The question of the real thrust characteristics 
of VRD can be answered by direct engine tests on test-beds 
to determine engine characteristics. Calculation of the actual 
combustion rules can have a substantial effect on the choice of 
engine type, parameters, and characteristics. In comprehensive 
analysis, calculation of the actual values of nr f° r an aircraft 
may provide results differing substantially from those of 
calculation analysis performed with npef = const. In particular, 
calculation of actual values may have a substantial effect on 
choice of operating modes (3) along the flight trajectory, 
choice of engine length, and thus choice of engine mass and 
amounts of combustor cooled surface, which is of fundamental 
importance with a limited cooling margin. When choosing opera- 
ting modes (3) it is also necessary to consider the possibility 
of providing (in these modes) cooling of powerplant and airframe 
components with hydrogen. 


The thrust of a PVRD depends on many factors, in particular: 

1 ( M £ ,fb « T.rgf » f m m u » . ^*o ) • ' 

The specific parameters are also functions of many f actors, /89 
especially for one fuel with a constant (invariable) flow duct: 


***** ~ s f ’ /*• T.r e f’ M c • ) 

s : J ( M f * T.ref 


Aircraft excess coefficients of fuel impulse and thrust depend 
on the aircraft's coefficient of external drag 
and are determined by relations 


Calculation of the actual values of hp e f allows choosing the 
correct cruise (Jr u = 0) and boost (Jr u > 0) modes. 

b) The Process in the PVRD and GPVRD Combustor 


Under conditions providing aerodynamic stabilization of a 
hydrogen flame or its self-ignition (at high Tg, corresponding 
to high flight M, and at adequate pressures in the combustor 
before fuel delivery), kinetic factors, especially a gas time- 
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of- stay in a combustor over 100 mm long, may not have a substan- 
tial effect on process efficiency. Even at gas velocities in 
the chamber W - 1500 m/sec with Z. R >.100 mm the ratio of the 
gas time-of-stay to the chemical reaction time (at local a = 1 
and T = T - 2500-2700) is much larger than unity and equals 


f - f 


15 - 20 


and the combustion efficiency is often substantially lower than 
1.0. We might assume that the combustion efficiency would be 
determined by fuel-air mixing. 


The efficiency of hydrogen's combustion with its diffusion 
mixing with air when delivered in sprays is determined basically 
by turbulent mixing processes that do not depend on the absolute 
values of the parameters in the chamber (Pr, W k ) or the chamber's 
absolute length, but are determined by the ratio of the cham- 
ber's length l ^ to the length of the diffusion flame in the 
duct l§, which (as in RPD [rocket-ramjet engines] [9]) is propor- 
tional to some complex > — the reduced relative length of the duct /90 
x 

Pr * „ P. Z 

1 - *pr • 

Analysis of the change in X pr from various relative parameters 
using turbulent- jet theory and also employing experimental regu- 
larities in flame-tongue change (in the accompanying flow of air 
and in unlimited space), obtained by Ferri-Zacchei [10] and 
corroborated by the author's experiments, revealed that for 
compressor less hydrogen VRD (PVRD, DGPVRD [dual-mode 
hypersonic ramjet engines] , GPVRD ) the reduced combustor length 
Xpr does not depend on the relative total outlet area of the 
fuel nozzles (F<j c /F k ) and does not depend on the relative tem- 
perature of the gaseous hydrogen fed into the chamber for com- 
bustion (0 T = Trj/T B ) . When providing equal "access" of air to 
each of nozzles in the fuel unit, if there is accompanying 
air (y = 0) or if it is fed into the chamber at an angle or per- 
pendicular to the engine's axis, the value of X pr is determined 
basically by the relation of the components (a or 3) and by 
the main geometric factor — the equivalent length of the 
combustor: 

fV7 , 

**eg 

where Dj( e q = /4 Fk/tt is the equivalent diameter of a combustor 
of arbitrary cross-section. 

The parameter that can be used to generalize the efficiency 
of the process in experiments with hydrogen PVRD and GPVRD with 
assorted designs is parameter X pr . The mathematical expression 
for the reduced length of the duct can be written as: 

*pr = * F (I> 
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where function F(a) is determined by the engine's operating con- 
ditions, bound into 3 groups in accordance with change in the 
flame plume length: 


I. < 1 
P UO * 





2.1 < oc < 2.5 * 3. at > 2.5 


The independence (for a hydrogen VRD) of the value of Xpr 
from the total area of the fuel nozzles (at Fr = const.) from 
air temperature (determined by the flight M) and fuel temperature 
(determined by its heating in the structural-component convec- 
tive cooling systems) is an extremely important circumstance 
permitting an arbitrary choice of nozzle area Ft c with variation 
of N c and experimental freedom from the requirement to provide 
simulation of Tg and T£. 

To obtain relation href = /(Xpr) and check the theoretical /91 
conclusions, a number of studies were performed using two-dimen- 
sional models (at y = 0) and with axially symmetric models. The 
studies were performed with external flow around the models of 
a hypersonic flow at Tft = 600-1500K and with an attached air duct 
at T* = 1000-1200K. 

Also employed in the analysis were foreign data from [4] 
and [5] (NASA and ONERA) at y = 90°, obtained for axially sym- 
metric PVRD with supersonic velocity in the chamber and mixed 
combustion (M R > 1) . 

Experimental analysis revealed the correctness of using 
Xpr as the modeling criterion. The obtained approximated (with 
10% precision) relation href = /(Xpr), shown in Fig. 3, is 
general and allows making a preliminary choice of combustor 
construction (N c , F^) , and also performing calculations of 

engine thrust and VKS boost characteristics considering the 
actual efficiency of engine operation. 


An approximation can be proposed for computer calculations: 
r) = 0.6 log Xpr. 

A number of studies (including those of NASA and ONERA) imply 
that increasing X above 25-30 with limitations Z- K /D e g = 2-3 
is very complicated because of the difficulty of cooling the 
combustor and nozzle, so the value X = 30 should be taken as 
a rough one when designing hydrogen PVRD (with subsonic, trans- 
sonic, and supersonic velocities in the combustors). 

c) Thrust Characteristics of VRD and Assessing VKS and TAKSi 
Boost Characteristics 


At X = 30 in PVRD and GPVRD there is combustion efficiency 
D ref = 0.95-0.65, depending on 3 (at 3 = 0.4-4), and regardless 
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of X at 3 = 1 there is a combustion efficiency "gap" of funda- /93 
mental importance for choosing engine operating modes (see below) ^ 
Figure 4 shows a comparison of calculated change in nref from 3 
with the experimental data of ONERA [4] , obtained in the M f = 5-7 
range. The calculated relation apparently agrees with the experi- 
mental data. This calculated relation is later taken for analysis 
of PVRD thrust characteristics (Fig. 4b, d) and aircraft excess /94 
thrust parameters (Fig. 4b) (with F . ^/F = 3) . 


When the aircraft is accelerated, to improve VRD efficiency 
with increasing speed (Mf) it is expedient to provide engine 
afterburning (increase 3); because of the combustion efficiency 
href "9 a P" at 3 = 1 and Mf - 6 it is expedient to shift "abruptly" 
to rich fuel-air mixture modes (3 > 1) • This peculiarity is 
apparently a general property of an actual engine, since it has 
already been noted when analyzing the acceleration of craft with 
an RPD . To provide for cruise flight (Jr u = 0) 
of a hypersonic aircraft at Mf = 6, after acceleration along 
line 3 opt at Jr Uit1? , v it is again necessary to throttle back 
"abruptly" from - 1.25 to 3 - 0.45 (see Fig. 4b). 

We shall examine the approximate possibility of accelerating 
a VKS with a PVRD. 


The PVRD fuel specific impulse and the fuel excess specific 
impulse Jr u in the aircraft system depends very highly on the 
flight M (Fig. 5) . In the Mf range from 2-3 to 10-12 in a hydro- 
gen PVRD at 3 = 1 = const, with 3 = 3 0 pt and J$u X ' the following 
relation can be employed: 


, rni nit .* 


■V 


( 2 ) 


where 5-2 and Jr u corresponds to the initial acceleration 
speed, which can be varied = 2-3) . 


Generally, with 3 = var. the value of 5 can be dependent 
on Mf; then in segments AM in the required range of the M interval 
we can make a numerical approximation 5 = /(Mf), assuming the 
end of the preceding segment as the beginning of the next. 

Without introducing a large error, we can assume 5=2= const, 
in the Mf < 12 range. This method allows substantially simplify- 
ing the estimate of the aircraft's acceleration characteristics 
in computations, as compared to 5 = var. 


From the general equation of aircraft motion, written as 

r " ~£r 

we can see that for rocket boost (when with a vertical launch 
sin v - 1, m/Ru << 1) and for aircraft boost (when m/Ru < 1, 
sin v - 0), for qualitative comparison under the same conditions 
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at- P*p*p4) 










the gravitational forces can be ignored. In this case, with /94 

conditions (2) for a craft with a VRD, considering that R u = 

GtJr u , it is easy to obtain a relation characterizing the craft' s /98 
accleration when burning on-board fuel. 


After integrating the equation of motion from its initial 
parameters to the current ones, we obtain: 


Y 


tV 


ini 




v init 


t» 


i 



(3) 


Relation (3) is general (for VRD and ZhRD) for the craft's 
boost segment where £ = const. In particular, for acceleration 
with a ZhRD (or an RDTT) , where £ = 0 (Jr u = const.), relation 
(3) is converted into the familiar equation of K. E. Tsiolkovskiy : 

* ' “init* *«. '• \ 

Equation (3), like Tsiolkovskiy ' s , is convenient in that, without 
construction analysis and without fairly complicated trajectory 
calculations, it allows estimating the main parameters, in par- 
ticular the required fuel reserve for accelerating a craft with 
a VRD. It is also convenient for comparative analysis of aero- 
space planes with ZhRD and VRD (PVRD) . 


A comparative evaluation of the acceleration characteristics 
of aircraft with ZhRD and VRD was performed in accordance with the 
characteristics of PVRD (Fig. 6) . The comparison revealed that 
with acceleration to Mf = 10 with a ZhRD at J Ru = 4000 N/kg-sec 
(JRu - 400 sec) the relative final mass of an aircraft m/mi n it 
should- ibe approximately 47%, and with acceleration of an air- 
craft with a TRD and hydrogen PVRD in the same velocity range, 
the final mass of the craft may be 78% of initial (22% of the 
initial mass of the craft should be fuel used for acceleration: 

5% kerosene for the TRD and 17% hydrogen for the PVRD) . 


Thus, we can state that the theoretical advantages of hydro- 
gen PVRD for boosting craft (of the VKS and TAKSi type) to 
M = 10-12 are indisputable. 
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NUMERICAL CALCULATION OF THE SEPARATION AND JOINING OF TWO- 
DIMENSIONAL SUPERSONIC FLOWS IN DUCTS WITH DISCONTINUOUS 

BOUNDARIES 

A. T. Berlyand and A. F. Chevagin 


The question of numerical calculation of a flow within a jet /100 
engine has been taken up in works dealing with development of 
such engines (GPVRD [hypersonic ramjet engines] , in particular) . 
Lectures devoted to elaborating the legacy and developing the 
ideas of F. A. Tsander have in past years reported on a series of 
works [1-3] dealing with study of the possibility of using for 
such calculations the previously suggested numerical method 
using automatic division of separations [4] . However, both the 
programs implementing this numerical method and other programs 
used for these purposes (see [5, 6, 7], for example) were 
generally intended for calculating the flow in a simple single 
duct. When developing GPVRD it often becomes necessary to 
examine a duct with different kinds of partitions, slots, branches, 
etc. Use of ordinary programs (of the type described in [6]) 
to calculate such flows makes every such calculation a unique 
and very laborious operation and does not allow carrying out 
parametric calculations in an acceptable period of time. In 
connection with the requirements for such calculations, it became 
necessary to create a multipurpose program implementing one of 
the familiar methods of calculating supersonic flows for regions 
with complex separation boundaries. 

This article deals with describing such a program imple- 
menting the numerical method in [4] and demonstrating its fit- /101 
ness using a number of examples. The choice of the calculation 
outline from [4J came about because it is specially oriented 
towards producing the most accurate solution near the solution 
features and, by omitting an outline of through calculation 
of the first order, it is devoid of a number of shortcomings 
present in ordinary through-calculation grid outlines [5] . 

An outline of the type in [4] can be used to construct gener- 
alized solutions of plane hyperbolic problems of gas dynamics 
in arbitrary regions for a fairly broad class of boundary con- 
ditions: solid boundaries, free surfaces, arbitrary incident 
flow, etc. [6] . The difficulties of using programs of the type 
in [5, 6] for calculating flows in separation boundaries is 
associated with the fact that we then lose the basic advantages 
of the through-calculation outline, the region under examination 
must be artificially divided into subregions, and the corres- 
ponding "joining" of solutions has to be performed, and this 
does not allow natural automation (an example of this kind of 
calculation is given in [6]) . It should be noted that for 
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sustained outlines of the type in [4] the automatic joining and 
division of flows should not encounter great technical (algo- 
rithmic) difficulties. When such outlines are employed, at each 
moment in time the program processes parameters of only one 
current section of the flow. 

Complete automation of sustained calculation of an arbitrary 
multiply connected region makes a very complex problem, generally 
leading to problems close in structure to those of image recog- 
nition theory. This apparently is because authors of the corres- 
ponding programs usually restrict themselves to calculating 
merely 1-connected regions and leave it to the user himself to 
analyze the flow in a complex region by employing the suggested 
program at some stages of this analysis (the exception is in 
periodic problems like the flow around grids of identical pro- 
files, which reduce to analyzing the flow within one period) . 

With this approach the user has to carry out both the part of the 
work actually difficult to formalize for computers, consisting of 
dividing the 1-connected subregions and of determining the se- 
quence of their calculation, which is usually not particularly 
difficult for man, and the purely technical part of the work, 
consisting of sorting the results of calculating each isolated 
subregion for use as initial data to calculate other subregions. 

The second part of the work, on the other hand, does not present 
any difficulties for the computer and is extremely cumbersome for /102 
man. This article offers a compromise approach to resolving this 
conflict. The user himself is assigned the division of a complex 
multiply connected region into a finite number of 1-connected sub- 
regions and determination of the order of calculating flow in 
these subregions. Figure 1 shows several examples of possible 
flow configurations which can be calculated using this program 
and means of breaking down the corresponding regions into sub- 
regions with the required numeration of their calculation sequences. 
Here the boundaries of all the subregions consist of solid walls 
and free surfaces. 

Information about the breakdown results is entered in the 
TABL secondary file, elements of which are entered into the 
computer at the very beginning of the calculations as part of the 
initial data. Using the TABL file parameters as an indicator, 
the program itself sorts the calculation results for each sub- 
region, and in those cases where they have to be initial data for 
calculation of other subregions it enters them into secondary file 
AOB, where they are stored until it becomes necessary to use them. 

For hyperbolic problems it is obvious that this breakdown is 
always possible, since for each subregion a mixed Cauchy problem 
is set up and there is strict demarcation and absence of reverse 
effects between the initial data at the left boundary of the sub- 
region and the calculation results at its right boundary. 
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KEY: (2) Initial data input 

(3) Generation of initial 
words of IOB subregion 

(4) Counting assembly 
[1-connected region] 


Fig. 2. 
(5) 
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(7) 


Sorting of counting 
results of subregion with 
number IOB 

(6) Sorting of counting 

results of subregion with 
number IOB 
Sorting of counting 
results of subregion 
with number IOB 
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A block diagram of the program is shown in Fig. 2. Its 
nucleus is the basic reference unit (block 4) , which consists of 
a program for numerical calculation of the flow in a 1-connected 
subregion with a sequential number (program identifier IOB) . In 
our case this block is practically equivalent to the program 
given in [6] . If necessary, it can easily be replaced by a sub- 
routine implementing any other calculation outline. The main 
program is made in the form of a "binding” of block 4, and into 
its function go preparation of initial data for solving the 
problem in the subregion with number IOB, triggering block 4, 
sorting of calculation results for this subregion, and entering 
of them into file AOB for storage and use in solution of the 
problem in subregions with numbers higher than AOB. 

Now a few words about the meaning and structure of the TABL 
and AOB files. Because of BESM-6 computer primary-storage limi- 
tations, the program is designed for solving problems in complex 
regions consisting of no more than 10 1-connected subregions. /105 
Files TABL and AOB are two-dimensional with the scale of TABL 
(1:10, 1:2) and AOB (1:10, 1:2100). Generally speaking, the 
values of TABL (IOB, 1) and TABL (IOB, 2) designate the number 
of the subregions in which the calculation results at the right 
boundary of the subregion with number IOB would be used as 
initial data. If TABL (IOB, 2) > 200, the value of TABL (IOB, 2) 
is not the number of a real subregion and is used as as indicator 
whose meaning will become clear from the following. The fact that 
in file TABL the second index of its elements assumes only values 
1 and 2 does not restrict the generality of the program and 
retains the possibility of calculating flow around grids con- 
sisting of several bodies, since we can always introduce an in- 
finitely small (not equal to computer zero) difference in X — 
the coordinates of the leading edge of each body forming the 
grid — and thereby introduce into the examination the corres- 
ponding number of additional theoretical subregions of infinitely 
short length (Fig. lg) . When two flows have to be united in 
any subregion, the calculation results for each of them should 
be transfered into the initial data for only one subregion in the 
upper or lower part of the corresponding initial layer. If the 
entire field of parameters in the last section of the subregion 
with number IOB is included in the field of initial parameters, 
in its lower part for the following subregion, the number of this 
subregion is entered by the user into TABL (IOB, 1) , and the num- 
ber 201 is entered into TABL (IOB, 2) ; if it is included in its 
upper part, however, the number of this subregion is also 
entered into TABL (IOB, 1) , but the number 202 is entered into 
TABL (IOB, 2) . Thus, in accordance with the values of the ele- 
ments of the TABL file, which by their meaning are only numbers 
of the corresponding subregions, during the calculation we can 
form elements of file AOB whose physical meaning is as profiles 
of gas-dynamic parameters at the entry to each of the isolated 
subregions. 
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When calculating subregions with number IOB the program 
operates in the following manner (see Fig. 2) . Before operation 
of basic computation block 4, block 3 is used to form a field of 
initial parameters at the input to the subregion with number IOB, 
which is taken as AOB [IOB, j] (j = 1-2100). After this, 
computation block 4 operates, which calculates the isolated 1- 
connected subregion. When the right boundary of the subregion 
with number IOB is reached, the calculation of this subregion is /106 
curtailed and the calculation results are sorted in accordance 
with the values of TABL [IOB, 1] and TABL [IOB, 2], In blocks 5 
and 7 the parameters obtained as a result of operation of block 4 
are sent to the appropriate place in file AOB, after which the 
number of the next calculated subregion IOB is increased by one 
(block 8) and the program sets about calculating the next 
subregion. 

With increase in the number of IOB and sequential calculation 
of the corresponding subregions, file AOB is gradually filled 
with the values of the initial parameters for calculating all the 
isolated subregions. On the other hand, as IOB increases, the 
values already entered into file AOB would gradually be used in 
calculating all the isolated subregions. This implies that the 
isolated subregions should be numbered by the program user so 
that all the subregions whose calculation results are needed for 
calculating subregions with number IOB would have a lower number 
than IOB. Then by the time that it is necessary to calculate 
the subregion with number IOB the corresponding elements of 
file AOB (IOB, j) would already be generated. Quite obviously, 
this kind of numeration is always possible. 

Let's examine numerical solution of some model problems. 

Figure 3 shows the solution of a problem whose structure is close 
to that of interference of flows between two profiles, given in 
[6], whose solution was obtained using "manual joining" of the 
solution for different subregions. In the interference problem 
examined in [6] it proved possible to calculate only the flow 
pattern between the profiles, because of the extremely cumbersome 
"manual" preparation of data for further flow calculation. 

The problem solution shown in Fig. 3 produced a detailed picture 
of the flow pattern in the region behind one of the profiles, 
the complexity of which is produced just by the mere presence of 
flow interference in the preceding subregions, each of which 
separately has a relatively simple structure. Similar calculations 
were also performed for a number of complicated configurations of 
solid surfaces, and the obtained flow patterns have a similar 
character. 

Another typical example of problems in which it is expedient 
to employ: this program is calculation of flow in a duct with 
perforated boundaries and a given pressure behind the perfora- 
tions. An example of calculating this kind of flow is shown in 
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Fig. 4. The characteristic feature of these problems is 
the fact that the boundaries of some of the isolated subregions /109 
partially are free surfaces. In distinction to the problems with 
free surfaces solved using the program in [ 6 J , in this case only 
a part of the surface is free and, as before, non-flow-through 
boundary conditions have to be set up for the remaining part. 

The program provided for special measures ensuring automatic 
transition at the necessary point from one boundary condition to 
another. Obviously, with this approach to calculation of flow 
in perforated boundaries the maximum number of perforation 
openings is limited to the maximum size of file TABL, since 
introduction of each new opening is associated with isolation 
of an additional subregion. It should be noted that problems 
concerning supersonic flow in perforated boundaries were quite 
successfully solved using the program given in [7] . There are, 
however, certain differences in the capabilities of these pro- 
grams. So in some cases it is convenient to use the program in 
[7] and in others to use this one. If there is a need for 
accuracy of calculations, characteristic of technical problems, 
and for evenness of incident flow, the program in [7] allows, 
without exiting to a peripheral storage device, calculating a 
flow with a large number of perforation openings. Because the 
program in [7] does not provide for use of conditions at a free 
surface, however, in those cases where a disturbance from the 
flow enters a perforation opening and may, by being reflected 
from the free boundary, affect the flow pattern within the main 
duct, a correct solution cannot be obtained using the program in 
[7] . A particularly unfavorable circumstance here is the fact 
that normally we cannot know beforehand whether or not such a 
situation arises. Our program is free of this shortcoming (see 
Fig. 4). Moreover, it can be used to calculate flow both above 
a perforation and below it, if all the flow remains supersonic. 

The convenience of employing this program for calculating 
two-dimensional double-flow air intakes of aircraft power plants 
is quite demonstrative. One possible design for such an air 
intake, widely used in the U.S.A. in designs for a hypersonic 
aircraft, is shown in Fig. 5. A power plant with the air intake 
shown in Fig. 5 provides acceptable characteristics in cruise 
flight at high supersonic and low hypersonic speeds (Mach 3-7) 

[8, 9]. 

Calculation of the flow in the lower duct, as in a normal 
air intake, can be performed both "manually ", by relations for 
shock waves in an inviscid flow, and based on many known numerical 
methods, and it obviously presents no additional difficulties. 

But calculation of the upper duct, which operates in an essen- 
tially uneven flow field behind a series of shock waves from 
the compression wedge which interact with rarefaction of waves 
from the partition separating the ducts, becomes complicated and 
cumbersome for both "manual" calculation and the familiar 
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numerical methods with "joining" of the solution at the boundary 
passing through the leading edge of the separating partition. In 
parametric calculations of flows in the ducts of this air 
intake, work would be extremely laborious and inefficient, 
depending on the air intake's geometry (compression wedge geometry, 
shape and location of separating partition and upper lip, etc.) 
in a wide range of flight Mach numbers and angles of attack. 

The complexity of flow calculation is significantly curbed 
using our program. Here it is sufficient to break down all of 
the flow field of interest to us into subregions and to generate 
file TABL in the following manner (Fig. 5): 


TML: 


3,2 

5,4 

2,202 

2,201 

4,201 

6,205 


The obtained calculation flow pattern in the air intake at 
M = 4 and angle of attack a = 0 is shown in Fig. 6 . The shock 
waves from the forward panel of the compression wedge pass in 
front of the edges of the lips of both air intake ducts, so 
the calculated flow factors through the air inlet ducts are 
fairly small: = 0.774, ^2 = 0.185. The flow in front of the 

inlet into the lower duct is fairly highly decelerated. The Mach 
numbers in this section vary from 1.832 to 2.683. Accordingly, 
the ratio of static pressure in the duct inlet section to the 
static pressure in the incident flow p/p varies from 4.925 to 
6.187. 

The upper duct operates in the flow field with local Mach 
numbers at the inlet exceeding the corresponding Mach numbers at 
the inlet to the lower duct, because of the ring of shock waves 
formed at the leading edge of the separating vane. For this /115 

reason, strong shock wave reflections from the lip form in the 
duct, which leads to high losses of total pressure. 

The distributions of pressure at the surfaces of both ducts 
are shown in Fig. 7 and 8 . Formation of shock waves at the lip 
(separating partition) in the lower duct is manifested in a sharp 
stepwise pressure increase in sections X = X^ and X = X 2 - A 
similar picture is also observed in the upper duct. The pressure 
on the panels of the compression wedge increases practically 
monotonically and in section X = X 3 (throat section) reaches the 
maximum value p/Poo = 21.5. Small changes in p/po, in sections 
X 4 and X 5 on the lip are evidence of the formation of weak shock 
waves . 
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ON THE THEORY OF AN ELECTROMAGNETIC RESONATOR ENGINE 


A. S. Dmitriyev 


In connection with the intensive development of laser and /117 
SVCh [microwave] technology, the literature widely discusses the 
idea of emloying a powerful source of electromagnetic radiation 
outside an accelerated craft for creating jet propulsion. For 
example, in [1, 2] there are examinations of reaction engine 
designs based on an evaporative mechanism of thrust, and [3] 
proposes accelerating a craft through the pressure of an expan- 
ding plasma obtained as a result of laser breakdown of the air 
and atmosphere near the craft. A diagram of a power plant com- 
bining the fundamental possibility of accerating an LA [aircraft] 
to relativistic velocities and a high factor of converting radia- 
tion energy into the kinetic energy of aircraft motion -- an 
electromagnetic resonator engine (EMRD) — is given in [4]. 

This article examines some aspects of this engine's theory. 

A schematic diagram of the EMRD is shown in Fig, 1. The LA 
is accelerated through the pressure of electromagnetic radiation 
on an open resonator, formed by mirrors 2 and 3 on the LA mirror 
(3) . The electromagnetic radiation is pumped by source 1 into 
the resonator through rectifier 4. Equilibrium pressure of 
electromagnetic radiation [ sia ] . Here 0 means the Q factor of 
the resonator. For an ideal system without loss, the Q factor 
is determined by the size of the resonator and the speed of the 
moving mirror. Acceleration of the craft continues asymptotically /119 
until the electromagnetic oscillations in the resonator are com- 
pletely damped after source 1 is switched off. When there is 
no side scattering and no losses in the mirrors and medium, 
the energy of electromagnetic oscillations is transformed into 
LA kinetic energy. During the transfer of energy the radiation 
entering the resonator gradually "reddens." 

The engine system provides for the presence in the fixed 
source and the LA of focused and strictly mutually oriented 
mirrors. This allows multiple use of the waves reflected alter- 
nately from each mirror to increase the LA momentum. It is 
through the multiple use of the momentum of photons imparting a 
small fraction of their energy to the LA with each reflection 
from its moving mirror that a high coefficient of converting 
electromagnetic oscillations into kinetic energy of the craft is 
achieved, wherein lies the critical advantage of EMRD over a 
photon engine. 
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One other advantage 
of EMRD is the high rate 
of LA acceleration at a 
given source power, as a 
result of the practically 
inertia-free conversion 
of one type of energy 
into another. 

For quantitative 
evaluation of EMRD effi- 
ciency we have to study 
the system's equation of 
motion. For the sake of 
simplicity, we shall 
examine a system without 
losses and introduce the 
following notations: 
m is the mass of the LA; 
X, F, and E are the size 
of the resonatorn the pressure force of radiation on 'the LA mir- 
ror, and the total energy of the electromagnetic field in the 
resonator, respectively. 

The acceleration of an LA is determined by relation: 

mi a F. ^ 



In turn, the pressure of radiation in the resonator is 
expressed by the total energy of electromagnetic radiation in the 
resonator and the size of the resonator 

F * E/x , (2) 

and since there is no loss in the system, the value of E at each 
point in time can be found from the law of conservation of energy 

E + \ mx* = ~ mil 4 E o 4 i Mt)At , (3) 


where Xq is the initial velocity of LA motion, Eq is the initial 
energy in the resonator; eg (t) is the power of the radiation 
source. Substituting E from (3) into (2) and expression (2) into 
(1) and introducing the notation 


cC * 4; mx* 4 E 0 

we obtain the desired equation of motion 
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mxii ♦«! mi ! 2 I (?) 

* t 

It is impossible to study analytically the nature of the 
solution of nonuniform nonlinear equation (4) . But analysis of 
the system's dynamics when the source is switched on, i.e., for 
equation 


mxii mx z - «c =0 (5) 

can be performed analytically using the phase plane method [5]. 
The solution of equation (5) is: 



where C is the integration constant. Consequently 



( 6 ) 


Thus, the phase path of equation (5) consists of two 
branches, and the actual segments of the branches would exist only 
at those values of X for which 


n n 


1 me 



> 0 


Considering that at the initial point in time X = Xq, we obtain 

C — _ * e 


The phase path degenerates at Eq > 0 into two straight lines 
parallel to the x axis. In this case the velocity of the craft 
is, of course, constant and independent of X (Fig. 2) . 

In the other physically realized cases Eq 0 [sic] and 
graphic representation reveals that the phase paths are open and 
for high X the velocity X asymptotically approaches the value 



It is interesting to note that in the absence of loss and 
with Eq > 0 in the system, both conversion of electromagnetic 
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radiation energy into kinetic 
energy of LA motion and the 
reverse process are possible. 
Here the kinetic energy of 
the craft is totally conver- 
ted into the energy of elec- 
tromagnetic radiation at 
that point of the phase 
path with coordinates 

* « x 0 /(JTLiL ), x *0. 

The nature of LA motion 
when the source is switched 
on with £g(t) = e 0 and t ^ ie 
dependence of the system's 
behavior on initial condi- 
tions can be obtained by nu- 
merically integrating equa- 
tion (4 ) . 

The behavior of the /122 
phase paths relative to the 
initial conditions is quali- 
tatively understood from Fig. 3. The closeness of the phase 
paths is evidence of the stability of the solutions with respect 
to the initial conditions. With the system parameters assigned 
for the graphs, in 15 min of acceleration the LA reaches orbital 
velocity. 

Of fundamental signi- 
ficance for the possibility 
if implementing EMRD is the 
law of radiation spectral- 
composition change with 
acceleration of the craft. 

We shall study the 
dependence of the radiation 
spectral composition in the 
resonator as a function of 
time, in the assumption 
of absence of interaction 
between the radiation of 
different frequencies. 

When a plane wave 
X km perpendicularly strikes 
a mirror surface moving 
at velocity, the frequency 


X km/sec 



ao *o it to 


Fig. 3. 
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changes in accordance with relation (the Doppler effect) : 


u* s. u) 


c- V 

C + V 


tjj (« 




(7) 


where w is the frequency of the incident wave; to' is the frequency 
of the reflected wave. One wall in the resonator is movable, 
and in time dt a wave with initial frequency oj strikes it 
(c dt)/2x times. Therefore, the frequency of an electromagnetic 
wave at time t + dt - u)(t + dt) is associated with the frequency 
of the wave at time t - co(t) by relation 

u (t 4 At) s ) 

or u)(t* At)-gnw(t) c » , „ 

ctt = 2* V c ** 


The left side of the last equation represents the total 
derivative of In f-u)(t)] with respect to time and, consequently, 
equals 




taJ'tt) 

"wTtT ' 


Thus, the equation for w ( t) is: 

. _S_ tn(\ - 1 )■ 18) 

u)(t) "2*. c 


Assume that at the initial point in time the spectral compo- 
sition of the radiation in the resonator is B(t) and no radiation 
arrives in the resonator from the source outside. For any fre- 
quency a) at the initial time, the law of its change over time is 
determined from equation (8) . At time t = 0 the spectral den- 
sity of B(co) corresponds to frequency to. We shall construct an 
equation for B(to[t]). Derivation of this equation is similar 
to that of equation (8). Upon reflection from the moving mirror, 
the amplitude of the original wave changes according to law /123 


E'*E(' 


1 - 


c 


i + -r 




C - V 
C + V 


-'l ~ *6 (- 


C - V 
C ♦ V 


)z 


19) 


Accordingly, the spectral density of the power, proportional 
to the square of the reflected wave's amplitude, equals 
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&' » 6 (1 - m ^ y. 


( 10 ) 


After transformations similar to those performed in the deri- 
vation of equation (8) , we obtain 


S‘(t) 

b(t) 



(II) 


Using the obtained equations we can calculate the spectral 
density of the radiation. To do this, we need to break down the 
spectral range into small segments and integrate equation (8) for 
the frequencies of all intervals with the initial conditions at 
the left end. We thereby find frequencies to . (t) corresponding 
to frequencies to^(0) at time t. 1 


w.(o) — u> t (t). 


Then, considering that the law of change of the square of 
the wave's amplitude (spectral density) with change in frequency 
is subject to equation (11) , we can construct a curve of the spec- 
tral composition at time t using the curve of spectral composi- 
tion at time t = 0. We can show that this process constitutes 
a solution using the method of characteristics (w(t) plays the 
role of characteristic) of some equation in partial derivatives 
of the first order for B(u>, t) . 


Note that along each curve lying in plane (t, co) we can 
treat co and B(co, t) as functions of t, so that the total deri- 
vative of function B(co, t) equals 


dft _ 9B 3 b dUo 

ett ' 9t + 3 u) ctt ' 


( 12 ) 


We shall now examine curve co (t) in plane (t, co) , which satis- 
fies equation (8). From (8) we can obtain explicit relation 
a)' (t) , which, for now, we shall conditionally write using 
co (t) , X, and X' 


u)'(t) 


8 2x 


tn (i -Z . 


(13) 


On the other side, on curve co ( t ) y we can obtain an explicit expres- 
sion for B' (t) from equation (11) . We shall write it as 

b'(t) = b(t) «n (<-«■•-!■ V (14) | 

Substituting (13) and (14) in (12) , we obtain the equation sought 
for B(co, t) 
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u)(tV c 
l • *(t) 




c 

2-x(t) 


tn(\ - 4 


X (t) V 
C ' 


(15) 


Considering that X/c << 1, equations (8) , (11) , and (15) can be 

simplified and written, respectively 


x(t> 

w ( t ) * x 1 1 > 


_ , x ( t ) 

b(t) * 4 x 1 1) 



w(t)x(t) 2x(t) 

-^TuT" s ’* (t) ”*7tT* 


(80 


CIIO 


(150 
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We examined the process of obtaining a relation when the 
source is switched on, i.e., the Cauchy problem for equation 
(15'). The problem with switching on of the source is a typical 
boundary value problem, for which 


8 ( w , t ) * »(w) at - vO > 0 , t %0 

) *■ ) at t > Ot 6) a u ) 9 


and the solution is sought in region w > 0, t > 0 [6]. Figure 5 
shows the diagram ^corresponding to (t, w) . The characteristics 
begin at positive semiaxis w and at positive semiaxis t. This 
solution can be linked with solution of the Cauchy problem in 
a manner based on the fact that this solution coincides with 
solution of the Cauchy problem, if the characteristic is con- 
tinued through points t = tg, w = wg to the w axis and the inter- 
section points are designated w = w(0). Then the boundary value 
problem is formulated as a Cauchy problem. 
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Analysis of the EMRD 
outline reveals that the 
basic parameters of the 
engine system are deter- 
mined by the characteris- 
tics of the mirrors, the 
radiation source, and the 
precision of stationary 
installation and LA orien- 
tation. In turn, the 
efficiency of the EMRD 
is determined most of all 
by the maximum distance 
of the craft — D — at 
which the coefficient of 
electromagnetic radiation 
energy conversion into 
LA kinetic energy is still 
adequately high and by 
the power of the power 
source. 

We can show that the maximum efficiency of energy conversion 
-- n -- between two mirrors through electromagnetic radiation 
depends only on parameter t: 

r 2 *3 / * , r 2 . 

where Rq and R 2 are mirror dimensions [7]. Relation n ( t ) for 
mirrors with round and square apertures is shown in Fig. 4. In 
the first case, lying along the x axis is value 

Z R, , 

and in the second t = (tt/4) (AD/L^L 2 ) • /126 

In both cases the parameter along the x axis can be written 
as x = 2 tt(XD/s), where s is the minimum total area of the receiver 
and transmitter mirrors' apertures needed to obtain the efficiency 
of transmission described by these curves. 

For x < 1.0 the transmission efficiency equals practically 

100 %. 


Q 



Fig. 5. 


With increasing distance, the efficiency of the EMRD falls 
in proportion to the reduction in resonator Q factor 
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The requirements for transmission efficiency are fairly strict. 
For example, with a 10% total efficiency of the system, the mini- 
mum allowable transmission efficiency is 99.9%. 

The prospects for creation of full-scale IEMRD should be 
linked first of all with work to create energy transmission 
systems. The combination of energy transmission systems with 
sail structures as a mirror on an LA may be the first full- 
scale realization of the idea of an electromagnetic resonator 
engine . 
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MODEL STUDIES ON THE DYNAMIC CHARACTERISTICS OF 
LIQUID FUEL ROCKET ENGINE PUMPS AND TURBINES 
IN TRANSIENT CONDITIONS 

V. M. Kalinin and V. A. Sherstyannikov 


Introduction. Launches of the first ZhRD [liquid fuel rocket/127 
engines] were worked out basically empirically, by testing various 
launch outlines and cyclograms directly in engine firing tests. 

This required large expenditures of time and resources. The ; 
existing launch calculation methods failed to account adequately 
for the basic factors affecting the process, and could not serve 
as a reliable tool for its development. Over 30% of the accidents 
with ZhRD took place during the launch phase. As more powerful 
and complicated engines were produced it became necessary to 
create more efficient and cheaper methods of trying out launches. 

This problem was sblved using the following procedures: 

1) studying the physical picture of engine start; 

2) creating methods of mathematical analysis of this process; 

3) developing methods of hydrodynamic modeling; 

Through development of the first and second areas, methods 
were mastered for high speed recording of basic engine parameters 
in start conditions, some important consistencies determining the 
chemophysical processes were discovered, and procedures for 
calculating launch using a computer were created which more 
fully accounted for the discovered consistencies. 

Work performed in the third area allowed creating methods 
of hydrodynamic modeling which provided for developing specific /128 
launch systems, especially for engines of complex design. Hydro- 
dynamic processes have a major role in ZhRD starting. This con- 
clusion was reached by almost all researchers investigating the 
operating process of actual engines. Hydrodynamic processes are 
virtually the only ones through which all other stages of the 
chemophysical process constituting ZhRD start can be affected. 

Thus, development of methods for studying and developing the ZhRD 
launch process traveled the road of thorough study of these pro- 
cesses and creating effective methods for their certain prediction 
and strict control. 

Development of hydrodynamic methods for modeling the launch 
process allowed significantly expanding the information about 
physical processes occurring during launch and creating an 
effective tool for prior working out of launch even before engine 
firing tests [1-3] . 
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The purpose of this article is to examine experimental 
methods based on hydrodynamic modeling of ZhRD systems and to 
describe the obtained results. 

Modeling Methods and Similarity Criteria 

Hydrodynamic modeling methods are used to study the space- 
time irregularity of engine accessory filling, irregularity of 
engine injector operation start, processes of emulsion feed of 
components, and the dynamic characteristics of pumps and turbines 
in start and transient conditions. 

The modeling consists essentially in the following. The 
assembly being studied (combustion chamber, turbo-pump assembly 
(TNA) , starting valve at the inlet to the engine, etc., for example) 
is mounted in the operating position in a hydraulic system with 
a modeling liquid, having fuel lines similar to the hydraulic 
lines of a rocket (Fig. 1) . The fuel tank pressurization and 
the resistance and inertia of the hydraulic system are estab- 
lished in accordance with the similarity criteria discussed 
below. The modeling liquid is delivered by switching on the 
starter and automatic start equipment through the control panel 
using a test bed allowing assignment of different times for 
shifting starter and launch device commands. All of the basic 
parameters needed for determining the dynamic characteristics of 
the assembly under study are recorded during testing using high /130 
speed measuring equipment. 



Fig. 1. 
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We use the following similarity criteria when modeling non- 
steady-state hydrodynamic processes: 

1) Complete geometric similarity of the hydraulic system 
(a full-scale assembly is employed) . 

2) Equality of piezometric head at the inlet to the assem- 
bly: 


( 






>«• 


3) equality of inertial and effective hydraulic resistance 
input lines: 


*• r <. 4 A • v * 


4) Equality of inherent self-oscillations of liquid in 
hydraulic system elements: 


( = (Wj , 

where Pj is the pressure of the liquid at the inlet, p £ is the 
density of the liquid, represents the lengths of the line 

segments, 3 j represents the areas of the passage segments of 
the lines, AP ^ j represents pressure losses in the lines, M is 
the model system, and H is the full-scale system. 

Meeting of these criteria provides identity of the temporal 
characteristics of the filling process, that is, the time needed 
for the components to reach the assembly and the amounts of 
fluid hammers and pressure dips in the full-scale and model sys- 
tems. When these conditions are observed, simple calculation 
using the density of the model liquid is used to determine the 
expected expenditures of actual components . 

Hydrodynamic tests are performed with thorough preparation 
of the assemblies being tested and detailed visualization of the 
liquid flow pattern (Fig. 2) . 

The firing assemblies, pumps, and TNA turbine are prepared 
with pressure, temperature, and phase sensors and transparent 
inserts for photographing the flow. Measured besides pressures 
and temperatures in the sections at the inlet and outlet of the 
assemblies are the pressures along the flow area (at the periphery 
of the spiral tap of the pump, in the radial clearance of the 
turbine, etc.). 

The combustion chamber cooling lines and fuel component 
feed lines are prepared with high speed sensors of the delivery 
of the component and small-calibrated pressure strain-gauge 
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Fig. 2. 


transducers. To visualize the process of operation start of 
the fuel injector, the combustion chamber is made with an in- 
verted injector — with external outlet of the liquid from the /132 
injector (Fig. 2) . This allows making motion picture films of 
the liquid emerging from the injector (in two projections) in 
the start process. 

Displacement of the controller plunger piston in the start 
mode is determined by inductive pickups mounted on its housing. 

To measure pressure drops and flows in feed lines for the fuel 
and oxidant of the combustion chamber and gas generator, differ- 
ential pressure pickups and induction turbine flowmeters are 
mounted . 

Hydrodynamic modeling methods have found extensive use. 

These methods are used to study the characteristics of non- 
steady-state processes and elaborate a number of assumptions for 
efficient organization of pre-firing processes. Some of the 
obtaiend results are examined below. 

II. Examples of Modeling 

1. Hammer during Pump Filling • > 

In the staring of engines in which the pumps are not filled 
with components before start, hammers accompanied by later wave 
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processes arise in the input hydraulic lines at the moment the 
pumps are filled. The allowable amounts of hammer are limited by 
the strength of the pump lines and housings. Modeling has revealed 
that two hammers develop when screw impeller pumps are filled. 

The first is associated with deceleration of the liquid in the 
screw, the second with generation of a critical flow condition 
in the volute diffuser. The second hammer is considerably greater 
than the first (Fig. 3) . Measurements have established that 
the basic part of a pump's hydraulic resistance in the filling 
mode is concentrated on the volute diffuser and amounts to approx- 
imately 85% of the total pump resistance. The nodal points on 
the resistance distribution characteristic along the pump tract 
(Fig. 3) correspond to: (1) resistance of the screw, (2) of the 

screw with impeller, (3) of the entire pump. Within each of 
these pump elements we can adopt a linear law of resistance and 
inertia distribution by volume. This relation can be used to 
calculate hammers for a whole series of pumps. The calculation 
results provide satisfactory agreement with experimental data. 
Because of this it is possible correctly to predict the amounts 
of hammer upon ZhRD start and adopt and approve in timely fashion 
measures for reducing them to a safe level. 
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2. Pressure and Torque of Pumps in Non-steady-state Conditions /134 

Non-steady-state processes upon start of ZhRD affect the 
operation of pumps basically as a consequence of the action of 
inertial masses of fluid in the flow part. In accordance with 
the equations of motion of the pump rotor liquid masses in its 
lines, the dynamic pressure and torque of the pump in transient 
modes are determined by expressions 



where H C ij. and M CT are the pressure and torque of the pump in 
the steady mode, -L*, . W, are coefficients characterizing 

the pump's inertia: 

v ■ | 

( l is the length of the path of movement of the liquid in the 
flow-through part of the pump, [illegible] is the area of the flow 
section of the element, z is the current value of the radius, 

6 is the angle between the vectors of relative and circumferential 
velocities, g is the acceleration of gravity, [illegible] is the 
liquid density, "K" is the impeller) . The above equations are 
approximate, since, of all the acting dynamic factors, only the 
inertia of the liquid and of the pump rotor is accounted for in 
them. 


Hydrodynamic modeling of pump operation in non-steady- state 
conditions has confirmed the need for employing these approximate 
relations to calculate the pressure and torque of pumps in ZhRD 
start modes. At relatively low rates of flow increase (dG/dt 
< 250 kg/sec/sec) the reduction of pressure in the TNA accelera- 
tion process is relatively low and does not exceed 5% of the pres- 
sure in the steady condition. At high rates, obtained with 
even opening of the pump valve (dG/dt = 2000-3000 kg/sec/sec) 
the drop in pressure reaches 30-40% (Fig. 4). Pressure drop 
occurs basically because of inertial losses in the diffuser. 

With increasing pump speed, the effect of liquid inertia on the 
pump decreases as a result basically of increase in the relative 
area of the diffuser throat. Acceleration of the rotor 
with real TNA acceleration rates has little effect on pump 
pressure. These conclusions have been tested on ZhRD pumps /136 

with specific speeds of n s = 36-75 at flow rate increases 
dG/dt < 3500 kg/sec/sec and rotor accelerations dn/dt < 10^ 
r .p.m./sec. 

The amount of torque on the pump shaft in transient modes 
(at dG/dt £ 500 kg/sec/sec and in the absence of cavitation) 
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Fig. 4. 


/135 


differ little from the amount of torque in the steady state /136 

(the difference is no more than 1%) . The obtained results lie 
at the foundation of mathematical modeling of ZhRD start. 

3. Dynamic Characteristics of Turbines 


Works studying the dynamics of ZhRD turbines in non-steady 
states have recently appeared in the literature. We have used 
the hydrodynamic modeling method to obtain some results related 
to i assessing the dynamic properties of turbines at high rates of 
change in gas flow and turbine rotation speed. The parameters 
of combined resistance of the circuit following the turbine were 
varied. Five different types of turbines in the TNA system were 
tested -- two centrifugal turbines with small tt t = 1.5-1. 8 and 
three partial axial-flow turbines with large tt-j = 15-30. Research 
revealed that in a wide range of values of parameter u/Cag the 
dynamic factors of rotating rotor acceleration have no marked 
effect on turbine efficiency. Thus, at accelerations dn/dt 
£ 2x10^ r.p.m./sec, characteristic of ZhRD start, the power char- 
acteristics of turbines can be considered quasi-steady. 

The volume aft of the turbine substantially affects the developed 
turbine power during ZhRD start and acceleration. Figure 5 
shows the consistencies in change of pressure in the gas duct 


121 




♦ 


OWGltfAL PAGE IS 
OF POOR QUALITY 



Fig. 5. 


from the turbine over time at different volumes aft of the turbine . /136 
For a sharp change in turbine operation conditions during testing, 
a fast-acting two-position valve was used which provides for 
increasing the input pressure with acceleration dPg/dt = (1000- 
1500) x 105 Pa/sec. Figure 5a shows the course of processes in 
a turbine with a short outlet line with exhaust to the atmosphere. 
Pressure was measured directly at the inlet to the turbine and 
aft of the impeller at its outlet. It can be seen that in this 
transient process, caused by an abrupt increase in pressure at 
the inlet to the turbine Pj*j, delay of pressure behind the turbine 
P 2 comprises a mere 0.002 sec. Similar results were obtained in 
tests of other turbines without a space aft of the turbine in a /138 
wide range of turbine operating-mode changes. The delay time here 
was 0.002-0.005 sec. 

Figure 5b shows the source of processes in tests of a turbine 
with a large volume aft of it in the form of a 60 L exhaust mani- 
fold at the outlet. In this case the delay time of the transient 
process for P 2 is ^0.45 sec, which substantially exceeds the 
values recorded in the absence of a volume aft of the turbine. 

Estimates were made of the effects of an exhaust manifold attached 
aft of the turbine on P 2 pressure transient processes. 

When calculating these processes we can use the continuity 
equation for the "turbine-gas duct" system with the assumption 
that this system consists of a simple pneumatic circuit comprising 
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the active hydraulic resistance of the turbine, a concentrated 
volume aft of the turbine, and a concentrated hydraulic resistance 
at the outlet from the volume. The results of calculations pro- 
vide satisfactory agreement with model processes (Fig. 5b) . This 
leads to the conclusion that in the low frequency region, char- 
acteristic of start processes and transient conditions of ZhRD 
(f = 0-100 Hz) , the size of the volume aft of the turbine is one 
of the basic dynamic factors determining the transient process 
in a system with turbine. 

4. Cavitation Interruptions of Pump Pressure 


Cavitation interruptions of pumps, caused by sharp pressure 
drops at the inlet to the pump, have a substantial effect on the 
process of starting and other transient conditions of power 
plants with ZhRD. In multiengine installations the engines are 
found to affect each other, a fact associated with hydraulic 
disturbances arising in the fuel supply system. When designing 
and developing starting systems we have to account for the 
dynamics of development of pump cavitation interruptions 
caused by induced vibrations, solitary pressure drops in inlet 
lines, or presence of gas inclusions in them. The hydraulic 
modeling method has allowed thorough study of the features of 
these highly unsteady processes. 

Figure 6 shows a typical picture of transient processes in 
a pump when a brief severe drop occurs in pressure at the inlet. 

With minor drops, the shape of the pressure curve behind the pump 

P 2 corresponds to the course of input pressure Pq, there is no 

phase shift, and the pump pressure does not change in the process 

of the transient condition. With severe drops the Pq and P 2 pres- /140 

sure change curves are different. Pressure curves for P K behind 

the impeller and P 2 behind the pump follow change in the curve 

of pressure Pq at the inlet with a delay and distortion of shape. 

Comparison of the nature of change in the integral flow 
rate difference 

and the reduced pressure head indicates that change in pump 
pressure corresponds approximately to the change in cavity volume. 
Results of tests on a wide range of pumps confirm the described 
consistencies and allow presenting a generalized dependence of 
pump pressure change on the volume of the forming cavity in the 
flow part of the pump H = f (v^) (Fig. 7). The points designated 
by O refer to a pump tested with delivery of gas inclusions to 
the inlet. Points marked x were obtained in tests with abrupt 
changes in resistance in the pump lines, and points marked X 
correspond to conditions of developed cavitation self —oscillations . 

The following approximations were obtained as a result of processing 
all the experimental data using the least-squares method: 

h - 1 + 0,0529 vl* - 0,8766 ** . 
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Based on the obtained generalized relation, mathematical /140 
models of the dynamics of cavitation interruptions of pump pres- 
sure were constructed: 


Model I 

Model II 

II 

a 2 

«4 

H , 


At 

f H ’ 



v : . 

-*1 

M 

j'(v R > ’ 

* 


j o 

Here ap h is the pump pressure according to the 

universal pressure characteristic; AP H is the pump pressure 
during cavitation; H = f (v K ) is the relative reduction of pressure 
during cavitation; v K = v/vg is the relative volume of the cavity; 
vg is the volume of the screw and impeller; G^ and G 2 represent 
the liquid flow at the inlet and outlet of the pump. The static 
cavitation (interruption) characteristics of pump pressure are 
used in the first model. 

r *t ( -rr£* ‘ “*§-)! and the variable of liquid stay time in 
the band t h = var. at G^ = G 2 • In the second model the static 
cavitation characteristics of maximum flow rates through the 

\ are used: 

6j at and V K »0 

“ & at & z > -6p r and = 0 

. 6 pr at > V„-»0 

The first model is valid for conditions with relatively minor /142 
cavitations (1 >. H _> 0.6), the second for the entire range 
(1 >_ H >_ 0) . Figure 6 gives the results of a comparison of ex- 
perimental transient processes with calculated ones. There are 
some differences in the quantitative indices of the processes. 
Introduction into the first model of correction SP^ = 0.1 x 10^ 

Pa for accuracy of measuring pressure Pp (t) provides better 
agreement. Thus, the obtained mathematical models allow accounting 
for the dynamics of cavitation interruptions of pump pressure in 
calculations of the start and other transient conditions of ZhRD. 

5. Diagnosing TNA Transient Conditions 

Technical diagnostics have been widely developed in aviation 
gas-turbine engines and have recently begun use in ZhRD. We used 
the following methods for diagnosing TNA in model conditions: 

— recording and analysis of TNA rotor mechanical vibrations 
produced by its hydrodynamic loads, mass imbalance, shaft 


pump 
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Fig. 8. 


elasticity and support pliability; /142 

— experimental calculation determination of hydrodynamic 
loadings of the rotor according to measurement of mean values 
of parameters of steady and transient conditions. To understand 
the complex interaction of hydrodynamic and mechanical forces, 
the pattern of combined action of the housing vibrations, radial 
and axial forces of shaft displacements, and the thermal state of 
the bearing are studied. The combination of measuring equipment 
implemented on a full-scale TNA includes strain gauges of radial 
forces, built into the bearing, a multi junction annular thermo- 
couple, current vortex sensors of shaft radial and axial displace- 
ments, vibration tranducers on the housing in the area of the 
support, and pressure-pulsation sensors at the pump's inlet and 
outlet. Figure 8 shows a diagnostic diagram of TNA preparation. 

Under modeling conditions, research is conducted by replacing 
the actual fuel components with water and the producer gas with 
air. Various methods of varying the loads on the TNA rotor are 
then used. A convenient means of simulating dynamic phenomena 
under these conditions is artificially bringing the pump into 
cavitation conditions. This method makes it possible, in 
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Fig. 9. 


/144 


particular, to change the hydrodynamic loads on the rotor within /142 
wide limits and to assign the desired level of modal parameters, 
including with respect to flow rate of cooled liquid through the 
bearing. This equipment allows obtaining complex information about 
the dynamic loadings and displacement of the rotor with errors /145 
not exceeding 10% in the 5-2000 Hz frequency range. 

Research has shown that significant increase in loads on the 
bearing and vibration displacements of the rotor are observed 
in transient conditions only with severe throttling of the pump 
— at low values of parameter Q/n (Fig. 9). Here even the direc- 
tion of the radial force is altered. If in condition Q/n nom. 
the force is directed towards large cross-sections of the diffuser 





then in throttling conditions there is a change in direction 
towards smaller diffuser cross-sections. At Q/n = 0.1 the 
change in force direction is ^50° . Cavitational self-oscillations 
in the "pump-inlet line" system do not lead to substantial 
additional loadings of the bearing. The obtained data are 
used to develop a mathematical model of a TNA as a diagnostic 
object. 
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HIERARCHY OF MATHEMATICAL MODELS OF 
AIRCRAFT CONTROL MOTORS 

V. N. Gladkova, V. V. Kokorin, Ye. V. Solov'yev 


Introduction. The modern approach to designing complex sys- /146 
terns, including motors (DU) for controlling aircraft (LA) , is 
based on comprehensive calculation of the interaction between 
both system components and external conditions. This leads to 
the need for discovering quantitative associations of 
DU operating processes with control system (SU) characteris- 
tics and the properties of the controlled motions of LA. This 
problem can be solved only using modeling methods, and particu- 
larly mathematical modeling methods, widely used in modern 
engineering practice. Mathematical modeling lies at the founda- 
tion of automated design, the methodology of which is based on 
use of a systems-technical approach. This approach, in turn, 
relies on the hierarchic principle, which provides a description 
of the operating processes of functional assemblies in complex 
systems with different degrees of simplification. Thus, for 
automated design of optimal DU for LA control we need to construct 
a hierarchic structure of a complex of mathematical models (MM) 
distinguished by a different level of simplification in describing 
DU function. This level is associated with the need for retain- 
ing in the MM those processes and phenomena selected as defini- 
tive during design. These processes have to be provided essen- 
tial similarity with the original and insubstantial difference 
of the model with respect to other indices of operating process 
quality which are not definitive. The limiting conditions in /147 
developing MM are cost and time factors. The former are deter- 
mined by the available computer facilities and the second by 
the time allotted to the design process. 

Problem statement and solution methods. The question of 
constructing a hierarchic complex of MM of the functioning of 
LA control system DU is associated most of all with providing a 
foundation for the choice of structure and parameters to iden- 
tify models ensuring that the modeled processes are equal to the 
redl ones. Assume that the real process of DU function and 
its model are characterized by change over time of value 
m(t) = {x 0 (t)>, where 0 = 1, 2, ..., L and a; M (t) = (x m g(t)}, 

where q = 1 , 2, ..., ..., I, which are mutually asso- 
ciated by relations F^[a:(t)] = Cj (t) , where j = 1, 2, ..., P 

and F MK [x M (t)] = C MK (t), where K = 1, 2, ..., P M , L and l are 

the total number of parameters characterizing the real process 
and its model, respectively, h is the number of parameters de- 
termining the existence of similarity between the studied pro- 
cess and the real ones; C . (t) and C.,„(t) are known parameters, 
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P u Pj 4 is the number of discovered interrelations between the real 
process and those processes implemented in the model, respectively 
(generally P^ << P) ; then the condition of adequacy of the modeled 
process to the real one can be written as: 

at * E, 


where 6^ isithe specific functional, 6$^ is the constraint, 
and e is a constant preassigned to the bounded quantity. 

The form of the specific functidnal and of the constraints is 
determined by the tasks of the modeled system in the design pro- 
cess and the characteristics of the computer support equipment. 
The form of the specific functional generally determines the form 
of the model and, in turn, the hierarchic structure of the 
complex of models is determined by the hierarchy of the 
specific functionals. 


When examining the DU as an element of the system controlling 
the motion of various LA relative to their centers of mass, in 
most cases it turns out that the operating process of aircraft 
control motors is pulsed in nature [1] . The pulsed nature of 
control motor operation is provided by SU command signals sent 
to the motor input with various durations t^ and various inter- 
vals t n (Fig. 1) . The sum t K + t n = T characterizes the period 
of alternation of pulsed motor actuations, and the inverse, 

1/T = f is the frequency of pulsed motor actuations. 



Fig. 1. 
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The output signals of the motor as an executive component /149 
in response to SU command are: the thrust momentum of the motor 


J . = /« P dt, characterizing the force of the motor on the LA 

go u 

and the: required fuel expenditure on its production 


T 

m = /- m(t) dt. In this connection, the idea of completeness 
of modeling the functioning of control motors has to be associa- 
ted with completeness of description of the thrust change over 
time — P(t) — and the fuel consumption m T (t). Since the 
actual output signals of the motor are shitted in time in 
models of DU functioning, we also have to account for lags of 
thrust — — t en — and cutoff t co (see Fig. 1) . The nature of 
change over time of P (t) , ■ At ( t) , and the lag times t en and t co 
are, to some extent, determined by the features of control motor 
construction and the type of fuel employed. Many works [2, 3] 
have shown, however, that the change of output signals over 
time is also largely determined by condition factors, i.e., the 
characteristics of combination of times tx and t n . When modeling 
this effect we have to examine both an account of the prehistory 
of the motor’s state, which considers the presence o<f heat and 
mass exchange between adjacent motor engagement pulses. This 
circumstance necessitates reproducing in the models the con- 
nection between the modeled values and the current values of 
tj^(t) and t n (t). Three characteristic cases determining the 
conditions of heat-mass exchange in pulsed motor operation con- 
ditions have to be examined; independent pulses (isolated 
engagements) , pulses associated by heat exchahge (shown in Fig. 1 
by solid lines) , and pulses associated by heat and mass exchange 
simultaneously (shown in Fig. 1 by solid and dot-and-dash 
lines) . 


The entire process of one act of control motor engagement 
can be represented as consisting of a minimum five consecutive 
time segments (Fig. 1): I — the segment of electromagnetic valve 
gear operation delay; II -- the segment of valve opening delay, 
fudl ignition lag, and non-steady thrust output to the nominal 
level; III — the segment of nominal thrust and fuel expenditure; 

IV — the segment of motor combustion chamber evacuation; V — 
the segment of fuel evaporation from the motor spaces aft of the 
valve . 

Depending on the duration of the input signals and their /150 
different combinations, in the process of DU functioning we 
may encounter different combinations of their linkage: I, II, III, 
IV, V, I, II ...; I, II, IV, I, II, ...; I, II, HI, I, II, ..., 
etc. Proper description of the changes over time of the modeled 
parameters P, At, t en , and t co in all possible combinations of 
input signals is provided by the structure of MM construction, 
a typical block diagram of which is shown in Fig. 2. The form 
of the approximating relations for each segment can be chosen in 
different ways. The degree of complexity of the description of 
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operating processes is determined by the required accuracy of /150 
reproduction of the modeled values and by the operating range of 
input signal combinations. The construction of MM of control 
motor function is based on the principles of structural and 
parametric identification of the modeled values to the real 
ones and is implemented using one method of regression analysis 
according to the outline shown in Fig. 3. The basic role in 
the process of identifying modeled and real processes is played 
by goodness of fit. 



Fig. 3. 

Analysis of the operating processes in control motors [2, 3] 
indicates that each engagement of a motor in the pulsed mode 
of operation for all possible combinations of command signals 
corresponds to a certain value of the specific thrust pulse, 
obtained in the form of ratio 

1 - k 

*imp‘ jJm(tVdLt 

Jyimp accounts for the force of the motor on the LA and the 
fuel expenditures needed to produce it. At the same time, the 
specific thrust pulse characterizes, with a given type of fuel, 
the degree of motor perfection. The specific thrust pulse ob- 
tained for any number of engagements (n) in a series of pulses 
reflects firstly the losses associated with the non-steady nature 
of triggering and cutoff processes, the proportion of which 
increases with diminishing duration of the command signal. The 
value Jyimp also accounts for heat and mass exchange between 
adjacent pulses, associated with the duration of the interval 
between engagements, that is, it accounts for the prehistory of 
the motor's state. Thus, Jyimp is determined by any of a set 
of mode factor combinations: Jyimp [t^/ t n , n] . The maximally 
possible value of the specific thrust pulse of a control motor 
is determined in the continuous mode of its operation (Jycont) • 
Its value is stable and determined by the type of fuel employed 
and the degree of mCtor construction perfection. The ratio 

J . /j ^ = J retains the properties of complex assessment 
yimp / ycont y 
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of the effect of the pulsed mode with respect to the ratio to /152 
the continuous mode and thereby characterizes the proportion of 
losses produced by the pulsed operating mode . Thus it is expe- 
dient to adopt a dimensionless complex as the goodness of fit of 
the mathematical model to a real item. 

The degree of completeness of the accounting in this cri- 
terion of the heat and mass exchange characteristics, associated 
with the prehistory of motor condition, can be used to con- 
struct its hierarchic structure, which reflects the different 
level of detail in description of operating processes in pulsed 
modes in mathematical models. 

Taken as lowest level of the hierarchy is the goodness of 
fit correponding to the condition of production . of independent 
pulses by the motor. Corresponding to this case is condition 
3y = F(tfc, n) ; for each single engagement n = 1. The MM of this 
level reproduce the continuous mode of operation and single mo- 
tor engagements between which there is no heat or mass exchange. 

The next higher level is taken as complication of the /153 

goodness of fit, which considers the’ connection between the 
specific thrust pulse and pulse mode parameters, Jy = F(t K , t n , n) 
under the condition that tjr = const., t n = const., and n n est* 

Thus, for each combination of t£ and t n the specific thrust pulse 
of the motor has a specific value with an established heat ex- 
change of the motor in the pulsed operating mode. The MM of 
this level are capable of reproducing established pulsed operating 
modes considering an established heat and mass exchange process. 

The next hierarchic level is associated with accounting for 
the goodness of fit of an unestablished pulsed mode or its 
non-steady state. An unestablished pulsed mode is manifested 
in the goodness of fit by accounting for the specific thrust 
pulse of the motor, associated with the sequential number of the 
engagement Jy = F(tK> t n , n) , under the condition that t^ = 
const., t n = const., and n >. 1. The non-steady nature of the 
pulsed mode is accounted for by the dependence of the goodness 
of fit on the values of t K and t n , which can be in any combina- 
tion at any point in time, but under the condition that each 
combination corresponds to a quasisteady pulsed mode. In this 
case Jy[tK/ t n , n] under the condition that tK = var.,t n = var. , 
and n n est . The MM of this level are capable of reproducing 
either unsteady pulsed modes, considering an unsteady process of 
heat and mass exchange, or an unsteady nature of change in 
established modes considering an established process of heat and 
mass exchange. 

At the highest hierarchic level in the goodness of fit 
we need simultaneous accounting of the effects of both the un- 
established pulsed mode and the unsteadiness of change in the 


135 



OWGfltAL PAGE IS 
GT POOR QUALITY 


(1) i 

UCrrtQVHU&t 

i 

uHCM&vauuu 


_ ( 5 ) 


isnjcruser- 
jt da»pmt 


ucntjrc hup 

r*ZKtTH&t 

j cS’pzynaS* 


ucntiranu* 

On&THbtX 

oSpasMoi 


4jc- 

aredotaTtfc 

cxut 

paabr* 



&£m** | 

2pe/rSar- IX -r/i.,/!) " i*ccowt ttc’Vat 

*eett/ i ./?*/ j /7»/7 est n ~ 1 

! z**Vaz , ±n*Vat ; H&J 


_U3) i (24) _ _ (25) ^gSjX^ifl (27) 

-- 

noun?**- j __ uMfiy/tbCHtw ycroHolutujutcfl U Hny/u>CHt>ix 

riuB \05 uhOwhoP perruM uMnyrtbCMaix pe-^u/^oS 

Hcae/ieu \ BKAtowenue j p**cunoS ( 26 ) 


ueoapxuvec- 

Wj aarf«| fif 

nadesru ~ 



Fig. 4. 


KEY: (1) Data source 

(2) Theoretical development 
and statistical data 

(3) Tests of mockups 

(4) Test of experimental 
models 

(5) Scientific research 

(6) Thrust characteristics 

(7) Calculated or experi- 
mental thrust 

(8) Integral thrust in 
established pulse mode 


(9) , (10) Instantaneous thrust 

(11) Consumption charateristics 

(12) Calculated experimental 
consumption 

(13) Integral consumption in 
established pulse mode 

(14) Integral consumption in 
different combinations 
of modes 

(15) Instantaneous thrust 
considering prehistory 

(Continued on next page) 


136 


















Caption for Fig. 4 (continued): 


(16) Types of mathematical models 

(17) simple 

(18) limited 

(19) incomplete a) 

(20) incomplete b) 

(21) complete 

(22) Goodness of fit 

(23) Areas of application 
for models 

(28) Hierarchic 


(24) Continuous mode. 

Single engagement 

(25) Established pulsed mode 

(26) a) Exit to established 
pulsed mode 

b) Combinations of 

established pulsed 
modes 

(27) Any combinations of 
pulsed modes 

level of model 


136b 



I 


unestablished inodes. In determining the criterion this appears /153 
as Jy[tK, t n/ n] , if t K = var., t n = var. , and n ^ 1. The MM of 
the highest level are capable of reproducing any combinations of 
pulsed operating inodes of a motor, considering heat and mass 
exchange processes. 

Results . The mathematiccal models are constructed in accor- 
dance with the hierarchy of the adopted criterion. Depending 
on the degree of completeness of ' description of pulsed mode 1 
characteristics in MM, models are called: simple, limited, incom- 
plete, and complete (Fig. 4) . 

In turn, a different degree of description detail of the /155 
modeled values is possible at each level of the hierarchy. As 
an example, Fig. 5 shows a comparative assessment of the lowest 
level MM with a different degree of detail in description of the 
modeled values with their reproduction of single engagements. 

In model 1 the processes of change in thrust and consumption over 
time are idealized by a relay function considering delay. In 
model 2 the process of thrust change is taken as the real one, but 
with a delay constant over time, and the change in consumption 
is idealized by a relay function. In Model 3 the process of 
thrust and consumption change is take as real with delay. In 
addition, a relation is introduced which reflects the connection 
between the process of pulse aftereffect with a prehistory of 
motor state with respect to the duration of the command signal. 
Experimental points obtained by calculation in [3] are given 
for comparison. We can see from comparison of the modeling results 
that the simplest models (1) can reproduce single engagements 
for a command signal duration of over 0.6 sec. At lower dura- 
tions the error of determining fuel consumption increases. Model 
2 has a systematic J^; reproduction error that increases as the 
duration of the command decreases , a fact showing the inadequacy 
of modeling only of real thrust. Model 3 indicates the modeled 
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Similar results are observed in comparative assessments of 
models of higher levels of MM hierarchic structure. /157 

Conclusion. The development of MM is associated first of all 
with choice of information sources. For the obtained hierarchic 
structure of MM the continuity is displayed both by the goodness 
of fit criteria and by the information sources (see Fig. 4). 

The goodness of fit criteria of each lower level are a special 
case of higher level MM criterion. The information used to con- 
struct the MM of higher levels includes the information sources 
of lower level MM. 

The cyclicity of MM development consists in the fact that 
information for higher level MM development can be partially 
obtained using models of the lower level. A shortage of informa- 
tion is filled in in the process of step-by-step DU design (Fig. 

6 ) . 



Fig. 6. 
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